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I.   SUMMARY  AND  PRACTICAL  INTERPRETATIONS 

The  results  of  this  research  did  not  completely  achieve  the  original 
objective  of  estimating  optimal  cropping  strategies  for  dryland  farming 
in  the  Northern  Great  Plains,  but  a  good  deal  of  progress  was  made  toward 
this  ambitious  objective.   Failures  at  various  locations  were  caused  by 
inadequacies  of  available  data;  usually  deficiencies  in  the  quality  of 
soil  moisture  measurements  in  the  experiments  were  the  problem.   But  during 
the  last  two  decades  there  has  been  an  apparent  reluctance  to  conduct 
controlled  experiments  over  an  extended  period  of  time  at  a  single  loca- 
tion.  The  research  reported  here  does  demonstrate  that  it  is  possible  to 
derive  economically  optimal  estimates  of  cropping  strategies  when  reasonably 
adequate  data  are  available,  experimental  data  which  is  very  feasible 
to  obtain  with  a  consistent  effort  directed  to  that  purpose. 

It  would  appear  that  quite  reliable  decision  rules  were  obtained  for 
the  Highwood  Bench-Fort  Benton  Area  of  Montana  from  only  seven  years  of 
experimental  data.   The  only  serious  weak  link  in  the  model  developed  for 
these  data  was  the  estimated  probability  distributions  for  plant  available 
soil  water  after  crops  of  winter  wheat  and  barley.   This  weakness  is 
unavoidable  with  only  seven  years  of  data  and  can  be  made  more  precise  by 
only  two  methods:   (1)  extension  of  the  experiment  over  additional  years 
which  is  more  costly  in  terms  of  the  delay  in  availability  of  the  informa- 
tion; and,  (2)  devote  a  considerable  research  effort  to  synthesis  of 
plant  available  soil  moisture  from  year  to  year  using  climatic  data  and 
evapotranspiration  formulae.   The  latter  is  no  simple  matter  if  a  useful 
level  of  accuracy  is  to  be  achieved. 


The  economic  analysis  on  the  Highwood  data  showed  that  the  alternate 
winter  wheat-fallow  rotation  is  an  excellent  cropping  system  if  saline  seep 
damages  to  the  land  are  not  a  problem  and  ij_  disease  and  weed  problems  can 
be  kept  under  control.   But  we  know  that  saline  seep  is  a  serious  problem 
in  this  region  and  it  would  be  irresponsible  to  ignore  land  losses  in 
order  to  receive  modestly  higher  annual  cash  flows. 

It  must  be  recognized  that  experimental  yields  do  not  always  reflect 
all  of  the  disease  and  weed  hazards  associated  with  various  crop  rota- 
tions because  sites  are  deliberately  chosen  to  avoid  such  problems  in 
order  to  prevent  the  confounding  of  these  problems  with  controlled  factors 
of  the  experiment.   Consequently,  we  would  expect  yields  from  the  alter- 
nate winter  wheat-fallow  rotation  in  the  experimental  data  to  be  overly 
optimistic  relative  to  rotations  which  included  both  fallow  and  barley 
with  winter  wheat  in  a  rotation.   Bringing  barley  into  the  rotation  would 
reduce  the  disease  and  weed  hazards  actually  faced  on  the  typical  wheat 
farm,  as  contrasted  to  the  selected  experimental  site.   In  view  of  this 
practical,  subjective  element  in  translating  experimental  results  to  on- 
farm  situations,  the  economic  analysis  for  the  Highwood  data  was  performed 
in  such  a  way  that  the  individual  farmer  can  make  the  subjective  evalua- 
tion for  his  particular  conditions,  possibly  with  advice  from  agronomists  or 
plant  pathologists. 

The  methodology  used  in  this  prime  winter  wheat  region  of  Montana  was 
to  estimate  the  expected  opportunity  loss  in  net  cash  returns  from  using 
various  probabilistic  strategies  which  would  tend  to  control  saline  seep. 
Then  these  expected  opportunity  losses  in  dollars  per  acre  were  crudely 


-3- 

cranslated  into  increased  yield  in  winter  wheat  production  required  to 
compensate  for  expected  opportunity  losses  for  a  given  strategy. 

One  of  the  most  promising  strategies  was  to  recrop  winter  wheat 
stubble  with  barley  in  the  spring  if  plant  available  soil  water  is  greater 
than  5.0  inches  and  to  recrop  barley  stubble  in  the  fall  with  winter  wheat 
if  plant  available  soil  water  exceeds  2.0  inches.   When  one  of  these  critical 
soil  water  levels  is  not  met,  the  land  is  fallowed,  and  winter  wheat  is 
always  planted  after  a  year  of  summer  fallow.   Average  land  use  frequencies 
in  this  strategy  are  approximately  25  percent  fallow,  46  percent  winter  wheat, 
and  29  percent  barley.   Fallow  is  chosen  in  this  strategy  only  when  the  soil 
is  quite  dry  at  the  time  a  decision  must  be  made,  which  gives  good  control 
of  deep  water  percolation,  at  least  on  the  average.   Control  of  deep  water 
percolation  provides  a  control  for  saline  seep. 

It  was  shown  that  less  than  a  bushel  per  acre  increase  in  yield  of 
winter  wheat,  which  can  be  attributed  to  getting  barley  into  the  cropping 
system  with  a  relative  frequency  of  about  30  percent,  will  justify  this 
cropping  strategy  over  alternate  winter  wheat-fallow.   However,  this  asser- 
tion does  not  take  into  consideration  differences  in  fixed  costs  under  the 
two  different  systems.   These  fixed  costs  and  any  other  problems  such  as 
labor  requirements  associated  with  the  change  from  winter  wheat-fallow  are 
quite  specific  to  individual  farms.   Some  farmers  may  already  have  slack 
in  their  machinery  and  labor  supplies  so  that  only  variable  costs  are 
relevant,  while  others  may  be  in  a  very  tight  situation. 

Using  the  budget  data  for  a  typical  farm  in  Chouteau  County  (Schaefer, 
et  al.  ,  1976c) ,  fixed  costs  for  machinery  would  be  about  $2  per  acre  more 
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under  the  above  described  strategy  than  under  alternate  wheat-fallow.   This 
measure  of  the  difference  was  obtained  by  calculating  the  weighted  average 
of  fixed  costs  in  machinery  for  land  in  fallow  and  in  crop  for  each  of 
the  two  cropping  systems;  the  weights  for  alternate  wheat-fallow  were  50-50, 
while  for  the  strategy  based  on  soil  water  they  were  25-75  for  fallow  and 
crop,  respectively.   This  approach  has  some  deficiencies  in  that  the  allo- 
cation of  fixed  costs  to  fallow  and  wheat  land  is  somewhat  arbitrary  in  a 
budget  study  such  as  in  (Schaefer,  et  al . ,  1976c),  but  the  estimate  would 
be  expected  to  be  reasonable. 

A  large  component  of  fixed  costs  is  in  harvesting  equipment.  It  would 
appear  to  be  economic  for  a  farmer  following  a  strategy  based  on  soil  water 
to  rely  partially  on  custom  harvesting  during  unusual  years  in  which  a  larger 
proportion  of  his  farm  was  in  crop.  A  strategy  which  gives  a  relative 
frequency  of  crop  equal  to  75  percent  will  involve  years  with  substantially 
more  or  less  than  75  percent  of  the  farm  in  crop  because  of  the  statistical 
nature  of  the  decision  rule  based  on  soil  water  at  seeding  time. 

Fluctuations  in  machinery  and  labor  requirements  can  be  substantially 
reduced  by  keeping  the  farm  divided  into  strips,  much  as  in  the  old  crop- 
fallow  system,  and  introducing  more  complicated  strategies  which  deliberately 
break  up  the  tendency  for  the  various  strips  to  be  in  crop  or  fallow  at 
the  same  time.   The  general  idea  is  to  use  different  critical  soil  water 
levels  in  the  decision  rules  for  two  separate  tracts  of  land  when  both 
tracts  are  in  the  same  land  use  (barley  or  winter  wheat  stubble).   This 
type  of  interrelated  strategy  for  the  two  tracts  will  put  one  tract  in 
fallow  and  the  other  in  crop  with  a  probability  equal  to  that  of  soil  water 


-5- 

falling  between  the  critical  soil  water  levels  used  on  the  separate 
tracts.   More  detail  on  these  general  strategies  is  given  in  (Burt  and 
Stauber)  . 

Plant  available  water  in  the  soil  is  not  often  measured  by  farmers 
and  would  entail  considerable  cost  if  a  good  sampling  procedure  were  used. 
As  a  practical  approximation,  depth  of  wet  soil  calibrated  to  the  soil 
texture  has  been  suggested.   Another  advantage  of  using  depth  of  wet  soil 
in  place  of  plant  available  water  is  that  this  substitution  can  be  used  as 
a  method  for  extrapolating  the  results  of  the  analysis  to  other  soils  than 
the  experimental  site. 

We  would  expect  the  qualitative  relationship  between  the  optimal 
critical  soil  water  level  in  a  decision  rule  and  the  water  holding  capacity 
of  the  soil  to  be  directly  related,  ceteris  partibus.   That  is  to  say, 
a  higher  soil  water  capacity  within  a  given  depth  range  would  imply  that 
the  best  choice  of  the  critical  level  for  soil  water,  above  which  a  crop 
is  planted  and  below  which  the  land  is  fallowed,  would  also  be  higher  and 
vice  versa.   By  a  higher  soil  water  capacity,  we  mean  relative  to  the 
experimental  site. 

In  view  of  this  direct  relationship,  a  proportional  adjustment 
between  the  critical  level  of  soil  water  and  field  capacity  of  the  soil 
would  seem  to  make  sense  as  an  approximate  adjustment  for  soils  not  too 
different  from  the  experimental  site.   If  depth  of  wet  soil  is  substituted 
for  plant  available  water,  calibrated  to  correspond  to  the  soil  at  the 
experimental  site  and  the  critical  level  of  a  decision  rule,  then  a  defi- 
nition of  the  strategy  in  units  of  wet  soil  at  seeding  time  will  automatically 
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make  the  proportionate  adjustment  wlien  the  strategy  is  apjiliod  to  other 
soils. 

Plant  available  soil  water  capacity  at  the  Highwood  Bench  experiment 
is  approximately  12.4  and  10.0  inches,  in  the  top  5  and  4  feet,  respectively. 
The  strategy  using  2.0  and  5.0  inches  of  water  as  the  critical  values  in 
the  fall  and  spring  would  transform  into  10  and  24  inches  of  wet  soil  in  the 
fall  and  spring,  respectively.   The  fall  decision  for  winter  wheat  considers 
soil  water  to  5  feet,  while  in  the  spring,  only  the  top  4  feet  of  soil 
water  is  relevant  for  barley. 

In  summary,  the  following  strategy  is  being  proposed  for  dryland  grain 
farming  in  the  Highwood-Fort  Benton  area  to  mitigate  saline  seep  problems 
with  a  minimum  of  financial  sacrifice: 

1)  Plant  winter  wheat  on  barley  stubble  if  depth  of  wet  soil 
measured  from  the  surface  is  greater  than  10  inches,  other- 
wise fallow  during  the  coming  year. 

2)  Plant  barley  on  winter  wheat  stubble  if  depth  of  wet  soil 
measured  from  the  surface  is  greater  than  24  inches,  other- 
wise fallow  during  the  summer. 

3)  Plant  winter  wheat  on  land  which  has  been  summer  fallowed 
regardless  of  soil  water. 

Average  annual  precipitation  at  the  experimental  site  is  about  15 
inches.   This  same  strategy  would  probably  give  good  results  in  slightly 
drier  areas  in  the  region  as  long  as  average  annual  precipitation  is  not 
less  than  14  inches  and  winter  wheat  is  still  adapted  to  the  climate.   How- 
ever, the  expected  opportunity  loss  would  be  higher  under  drier  climates 
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aiid  the  strategy  should  maybe  be  adjusted  to  somewhat  greater  depths  of 
wet  soil.   For  example  in  areas  of  only  13  inches  average  annual  precipi- 
tation the  critical  levels  10  and  24  inches  could  logically  be  replaced 
by  say  14  and  28  inches  of  wet  soil,  but  we  have  no  objective  basis  for 
such  adjustments. 

Results  for  the  Spring  Wheat  Region  of  Northern  Montana  and  North 
Dakota  suggested  that  too  much  fallow  is  being  practiced  regardless  of  any 
saline  seep  problems  which  might  be  present.  Under  a  simple  long-run 
average  net  returns  criterion,  the  results  suggested  that  the  very  eastern 
part  of  Montana  and  all  of  North  Dakota  north  of  Mandan  should  probably  be 
in  continuous  cropping,  without  any  fallow  except  for  cases  where  benefits 
from  weed  control  can  justify  occasional  summer  fallow. 

There  was  a  large  data  gap  between  Eastern  Montana  and  Havre  in  North 
Central  Montana,  so  we  have  a  large  expanse  of  spring  wheat  country  where 
nothing  definitive  can  be  said.   The  optimal  strategy  at  Havre,  using  long- 
run  average  net  returns  in  a  simple  model,  was  to  plant  wheat  if  plant 
available  soil  water  is  greater  than  about  3.5  inches  and  fallow  otherwise. 
The  critical  level  of  soil  water  was  between  3.0  and  4.0  inches  depending 
on  wheat  price  and  yield  assumptions.   Data  did  not  permit  an  analysis 
of  other  crops,  such  as  barley,  in  the  strategy. 

The  economic  analysis  showed  that  using  3.0  inches  as  the  critical 
value  resulted  in  a  very  small  expected  opportunity  loss  and  this 
strategy  appears  to  be  quite  effective  for  saline  seep  control.   First 
of  all,  the  land  is  in  crop  about  two-thirds  of  the  time  under  this  strategy, 
But  just  as  important,  the  land  is  fallowed  only  when  soil  water  at  seeding 
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time  in  the  spring  is  less  than  3.0  inches.   By  contrast,  the  land  is 
in  fallow  half  of  the  time  when  alternate  wheat- fallow  is  practiced,  and 
during  the  years  of  fallow,  soil  water  in  the  spring  is  greater  than  3.0 
inches  64  percent  of  the  time.   Both  of  these  considerations  make  the 
strategy  based  on  soil  water  much  better  from  the  pciiU  of  view  of  saline 
seep  control. 

Transforming  the  strategy  to  inches  of  wet  soil  measured  from  the 
surface  gives  a  critical  level  of  18  Inches,  which  we  propose  as  a  practical 
strategy  for  the  spring  wheat  region  around  Havre,  Montana.   Average 
annual  precipitation  for  the  Havre  site  is  about  11.6  inches  and  average 
annual  precipitation  at  Sidney,  Montana  where  continuous  cropping  was  optimal 
is  about  14  inches.   It  would  appear  that  fallow  should  start  coming  into 
the  optimal  strategy  somewhere  between  11.6  and  14  inches  of  average 
annual  precipitation,  but  we  are  working  with  quite  limited  information 
and  average  annual  precipitation,  even  at  the  same  latitude,  does  not 
tell  the  entire  story  on  climate. 

Nevertheless,  let  us  make  two  simple  assumptions  and  estimate  a 
"reasonable"  strategy  for  the  spring  wheat  region  of  Montana  where  the 
critical  value  for  inches  of  wet  soil  at  seeding  time  is  correlated  with 
average  annual  precipitation.   First,  we  assume  that  fallow  should  start 
to  enter  the  optimal  strategy  at  about  14.0  inches  average  annual  precipi- 
tation.  Second,  we  assume  that  the  critical  value  for  inches  of  wet  soil 
changes  linearly  with  average  annual  precipitation. 

These  two  assumptions  give  us  a  linear  relationship  between  two 
variables  and  two  points  on  that  linear  relationship,  one  point  being 


associated  with  Havre  and  the  other  determined  by  our  first  assumption. 
The  linear  relationship  implied  by  our  assumptions  is  presented  in  tabular 
form  in  Table  A.   We  even  extrapolated  to  lower  precipitation  levels  than 
Havre  at  the  top  of  the  table.   Recall  that  the  critical  level  for  soil 
water  at  Havre  was  adjusted  downward  somewhat  from  the  "optimal"  level 
based  on  average  annual  net  returns  in  order  to  improve  saline  seep  con- 
trol.  Therefore,  the  strategies  in  Table  A  are  relatively  "dry"  ones  to 
help  control  saline  seep. 

A  more  "conservaL  i  vc"  sLrategy  In  that  it  dc-vi^ites  less  from  the 
traditional  rotation  of  alternate  crop-Laliow  would  be  to  use  4.0  inches 
of  plant  available  water  as  the  critical  value  at  Havre.   If  this  is  done 
and  the  critical  level  transformed  into  inches  of  wet  soil  measured 
from  the  surface,  the  counterpart  of  Table  A  is  Table  B.   If  no  saline 
seep  problems  are  being  experienced,  or  salts  contamination  of  ground  water 
supplies  is  no  consideration,  then  the  strategies  implied  by  Table  B 
are  superior  on  economic  grounds. 

The  critical  levels  given  in  Tables  A  and  B  are  potential  approxima- 
tions for  barley  as  well  as  spring  wheat,  although  there  are  obvious  risks 
in  this  kind  of  extrapolation.   Since  barley  does  not  root  quite  as  deeply 
as  spring  wheat,  the  critical  levels  might  be  adjusted  downward  by  about  10 
to  20  percent.   One  result  of  this  study  which  casts  considerable  doubt  on 
using  Tables  A  and  B  for  barley  is  the  "pure  fallow"  effect  on  yield  obtained 
for  barley  planted  on  fallow  in  the  Highwood  Bench  experiment  (see  Section  II) 
If  this  "pure  fallow"  effect  is  greater  for  barley  than  spring  wheat,  then 
a  reduction  in  the  critical  levels  of  Tables  A  and  B  would  be  in  the  wrong 
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T;ible  A 


Critical  Depths  of  Wet  Soil  for 
Spring  Wheat  Strategies 


Average  Annual  Critical  Depth 

Precipitation  of  Wet  Soil 

(inches)  (inches) 


11.00  22.5 

11.25  20.6 

11.50  18.8 

11.75  16.9 

12.00  15.0 

12.25  13.3 

12.50  11.3 

12.75  9.4 

13.00  7.5 

13.25  5.6 

13.50  3.8 

13.75  1.9 

14.00  0.0 


-11- 


Table  B 


Conservative  Critical  Depths  of  Wet  Soil 
For  Spring  Wheat  Strategies 


Average  Annual  Critical  Depth 

Precipitation  of  Wet  Soil 

(inches)  (inches) 


11.00  30.0 

11.25  27.5 

11.50  25.0 

11.75  22.5 

12.00  20.0 

12.25  17.5 

12.50  15.0 

12.75  12.5 

13.00  10.0 

13.25  7.5 

13.50  5.0 

13.75  2.5 

14.00  0.0 
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direction.   Obviously,  we  are  not  in  a  position  to  say  much  about  an 
optimal  strategy  for  recropping  with  barley,  more  experimental  data  is 
needed. 

One  additional  warning  is  in  order  for  both  the  Spring  Wheat  Region 
and  the  winter  wheat-barley  strategies  for  North  Central  Montana.   In  both 
cases,  empirical  information  was  from  open  soils  where  there  were  no  root 
restricting  layers.   Any  root  restricting  layers  will  reduce  the  soil 
water  reservoir  available  for  storing  water  during  a  year  of  fallow.   The 
practical  versions  of  the  strategies  where  the  critical  levels  are  expressed 
in  inches  of  wet  soil  will  partially  adjust  for  this  possibility,  but  the 
adjustments  are  much  more  appropriate  for  variations  due  to  soil  texture 
than  for  root  restricting  layers. 

Most  of  the  benefits  of  fallow  emanate  from  increased  water  storage 
in  the  deeper  layers  of  the  soil,  say  between  two  and  four  feet  from  the 
surface,  except  in  extremely  dry  years.   If  these  lower  depths  for  storage 
are  unavailable,  fallow  would  be  beneficial  only  in  exceptionally  dry 
years.   A  good  example  of  this  condition  is  much  of  the  Judith  Basin  in 
Central  Montana  where  a  gravelly  layer  restricts  root  penetration  at  about 
24  inches.   Experimental  results  show  hardly  any  benefits  from  fallow  wliich 
can  be  attributed  to  soil  water  storage  in  this  area,  but  average  annual 
precipitation  exceeds  15  inches  in  most  of  the  area  too.   On  the  basis  of 
results  from  this  study,  we  would  expect  fallow  to  be  of  questionable  value 
for  water  storage  in  such  a  soil  situation  even  with  annual  average  precipi- 
tation as  low  as  13  inches. 

A  comparison  of  two  maps  of  Montana,  one  showing  saline  seep  problem 
areas  and  the  other  average  annual  precipitation  contours,  is  quite  suggestive 
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that  13  to  14  inches  average  annual  precipitation  is  a  necessary  condition 
for  saline  seep  problems.   The  precipitation  map  had  contours  by  two-inch 
intervals  so  that  the  comparison  had  to  be  rather  crude,  but  a  region  of 
average  annual  precipitation  with  14  or  more  inches  can  be  identified  for 
most  of  the  major  saline  seep  areas  in  the  state. 

The  cropping  strategy  analysis  suggested  that  on  purely  economic 
grounds,  fallow  should  be  practiced  to  a  limited  extent  in  the  Spring  Wheat 
Region  when  average  annual  precipitation  exceeds  13  inches  (see  Tables  A 
and  B).   Consequently,  the  saline  seep  problems  in  the  Spring  Wheat  Region 
would  seem  to  be  largely  solved  by  merely  maximizing  average  profits  with- 
out regard  to  the  saline  seep  problems;  clearly  a  case  of  no  sacrifice  to 
deal  with  saline  seep.   Fallow  would  be  used  only  as  needed  to  deal  with 
weed  problems  or  an  occasional  extra  dry  soil  at  planting  time. 


II.   THE  HIGHWOOD-FORT  BENTON  REGION  OF  MONTANA 

This  is  one  of  the  most  productive  agricultural  areas  of  Montana, 
and  also  a  major  saline  seep  probJem  area.   Climate  and  soils  are  excel- 
lent for  winter  wheat  and  other  cereal  grains.   The  predominant  rota- 
tion since  World  War  II  has  been  alternate  winter  wheat-fallow.   Consider- 
able research  effort  has  been  devoted  in  recent  years  to  analyze  the 
saline  seep  problem.   The  concensus  of  opinion  is  that  the  primary  source 
of  water  flow  for  the  saline  seeps  emanates  from  deep  percolation  of 
water  through  the  crop  root  zone  during  the  summer  fallow  period. 

The  obvious  solution  to  the  problem  is  reduced  summer  fallow  in 
the  cropping  systems  followed  in  the  area.   The  research  being  reported 
here  was  an  effort  to  determine  the  economic  incentives  for  continued 
summer  fallow,  and  analyze  alternative  cropping  strategies  which  would 
reduce  deep  water  percolation. 

The  only  useful  empirical  data  available  is  a  seven-year  experiment 
which  has  been  performed  on  the  Highwood  bench  (1970-76).   Tlie  experi- 
ment was  organized  and  started  in  1969  by  Ray  Choriki,  and  in  recent 
years,  has  been  conducted  by  Paul  Brown  with  the  Agricultural  Research 
Service,  USDA. 

The  data  analyzed  were  from  three  rotations: 

1.  Barley-Winter  Wheat 

2.  Fallow-Barley 

3.  Fallow-Winter  Wheat 
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Each  of  these  rotations  was  used  in  an  experiment  with  nitrogen  fer- 
tilizer as  a  controlled  factor.   Data  from  four  levels  of  nitrogen  were 
used,  0,  30,  60,  and  90  pounds. 

Plant  available  soil  water  was  measured  to  five  feet  in  the  fall 
for  purposes  of  winter  wheat,  and  to  four  feet  in  the  spring  for  use 
with  barley.   The  silty  clay  loam  has  a  relatively  large  storage 
capacity,  about  12.4  and  10.0  inches  of  plant  available  water  to  five 
and  four  feet,  respectively.   Average  annual  precipitation  for  the  area 
is  a  little  under  15  inches. 

The  Statistical  Model 
A  yield  response  surface  for  each  winter  wheat  and  barley  was 
estimated,  where  the  primary  variables  were  total  water  available  to 
the  crop  and  nitrogen  fertilizer.   Total  water  was  approximated  by 
plant  available  soil  water  at  seeding  plus  the  sum  of  precipitation 
during  May,  June,  and  July.   Nitrogen  applied  as  fertilizer  was  used 
directly  when  the  crop  was  planted  on  stubble,  but  when  planted  on 
fallow,  a  parameter  was  added  to  applied  fertilizer  nitrogen  to  account 
for  extra  nitrogen  mineralized  in  the  soil  during  the  summer  fallow 
year.   This  parameter  was  estimated  statistically  along  with  other  para- 
meters of  the  model,  and  is  indirectly  determined  by  yield  response 
across  nitrogen  levels  of  the  experimental  data. 
The  following  notation  is  introduced: 

N  =  total  nitrogen  measured  on  a  scale  with  origin  (measure 

zero)  at  that  level  associated  with  no  applied  fertilizer 
for  a  crop  planted  on  stubble  (pounds  per  acre) 
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W  =  plant  available  soil  water  at  seeding  plus  precipitation 

during  M.iy  through  July  (inches) 

fO    if  crop  is  planted  on  stubble 

[JL   if  crop  is  planted  on  fallow 

y  =  crop  yield  (bushels  per  acre) 

The  postulated  functional  form  for  the  yield  response  surface  is 

(2.1)   y  =  b^  +  b,N  +  b^N^  +  b-N^  +  b.W  +  b^W^  +  b^W"^  +  b^NW  +  b„N^W 
01/      J      H  5  6      7      8 

+  b  NW^  +  b  X, 
9       10 

N  =  (applied  nitrogen)  +  a 

The  parameters  to  be  estimated  statistically  are  b„,  b,  ,  .  .  .  b-,^  and 

a.   The  last  parameter,  a,  is  interpreted  as  the  amount  of  nitrogen 

mineralized  during  the  fallow  year.   The  parameter  biQ  is  a  constant  which 

is  added  to  yield  when  the  crop  is  planted  on  fallow  (see  the  definition 

of  X).   The  remainder  of  the  yield  equation  is  a  general  polynomial  in 

9  o 

N  and  W  of  third  degree.  The  terms  in  NW,  N  W,  and  NW  are  interactions 
which  jointly  reflect  the  nonadditivity  of  the  net  effects  from  nitrogen 
and  water.  The  polynomial  form  for  the  response  surface  is  justified  by 
Taylor's  formula  of  elementary  calculus;  any  analytic  function  can  be 
approximated  by  a  truncated  polynomial  containing  a  sufficient  number  of 
terms. 

Let  the  response  function  in  (2.1)  be  written  in  abbreviated  form 
as  f(N,  W,  X).   Then  the  statistical  specifications  used  to  fit  the 
equation  to  the  data  are 


-17- 


(2.2)   y   =  f(N   ,  W   ,  X   )v.       ,     j  =  1,  2,  ...  8 
Jt^      Jt^   ^^        J"^   J^  t  =  1,  2,  ...  7 

The  subscript  t  is  associated  with  a  year  of  the  experiment,  and  j  desig- 
nates a  treatment  combination,  namely,  a  level  of  nitrogen  fertilizer 

applied  and  the  crop  planted  on  either  stubble  or  fallow.   Thus,  y    is 

jt 

observed  crop  yield,  while  N   ,  W   ,  and  X   are  the  associated  indepen- 

jt    Jt       Jt 

dent  variables  for  a  given  observation. 

The  multiplicative  variable  v.   is  the  statistical  component  of  the 
model,  an  error  term  if  you  please.   There  are  many  sources  and  dimensions 
to  this  error  term,  such  as  measurement  error  in  harvesting  the  plots, 
soil  variation  over  plots,  plant  available  water  variation  over  plots, 
variation  in  stand,  etc.   All  of  these  sources  of  error  are  random  across 
treatments  for  a  given  year's  data  because  of  the  randomized  experimental 
design. 

However,  another  source  of  error  which  is  year-specific  is  variation 
in  climatic  conditions  from  year  to  year,  other  than  the  variation  captured 
by  the  water  variable  W  included  in  the  equation.   These  other  sources  of 
variation  are  such  things  as  temperature  experiencad  before  and  during  the 
growing  season,  the  distribution  of  precipitation  during  the  season,  winter 
and  early  spring  precipitation  before  May  1,  disease  and  insect  damage, 
relative  humidity  and  wind,  and  probably  many  other  things.   The  important 
aspect  to  recognize  with  regard  to  the  year-specific  component  of  the  error 
term  is  that  all  treatments  in  a  given  year  experience  the  same  random 
shock. 

In  view  of  the  above  discussion,  it  makes  sense  to  think  of  the  error 
term  as  comprised  of  two  factors. 
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(2.3)  V.   =  e,    u,, 

where  e.   is  a  purely  randim  factor  associated  with  normal  exper imcMiLa 1 
unit  variability  In  a  KLVcn  year,  and  u  ^^  is  a  random  shock  imposed  on 
yield  from  the  many  infiuences  which  are  specific  to  a  j^iven  year. 

The  multiplicative  error  term  is  justified  on  agronomic  theory 
considerations  as  well  as  empirical  results;  the  alternative  would  be 
an  additive  error. 

The  model  of  (2.2)  is  transformed  into  an  additive  error  model  by 
taking  the  natural  logarithm  of  both  sides  of  the   equation  to  get 

(2.4)  log  y_.^  =  log  (f(Nj^,  W^^,  X.^))    +   log  v.^ 

Using  (2.3), 

(2.5)  log  V   =  log  E.   +  log  u^ 

Letting  log  v   be  denoted  by  v'  ,  and  recalling  the  above  discussion 
Jt  Jt' 

of  the  two  error  components  l.   and  u    the  statistical  specifications  for 

^jt  ^"  = 

Cov(v'  ,  v'  )  =  a 

(2.6)  ^^        J^     ^J 

Cov(v;  ,  v'.  )  =  0   ,   t  5^  m 
It    jm 

We  actually  make  a  somewhat  stronger  assumption  than  (2.6)  based 

on  the  following  argument.   Since  we  are  dealing  with  data  from  a  randomized 

experimental  design,  the  following  relationship  holds  in  the  correlations: 


(2.7)   f... 
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In  addition,  we  assume  that  the  standard  deviation  of  y.   is  proportional 

to  the  mean  of  y^j-,  which  implies  that 

(2.8)   Oii  =  o^  ,    i  =  1,  2,  ...  8 

Assumptions  (2.7)  and  (2.8)  together  imply  that  there  is  only  one  extra 

unknown  parameter  in  tiiis  model  comijared  to  the  error  specification  in 

the  classical  regression  model,  i.e.,  the  correlation  parameter  for 

treatments  in  a  given  year,  designated  as  p  in  (2.7)  above. 

There  is  a  complication  in  estimating  the  parameters  in  (2.4)  since 
the  parameters  enter  nonlinearly  in  a  least  squares  criterion.   This 
problem  is  handled  by  a  generalized  nonlinear  least  squares  algorithm 
developed  by  one  of  the  authors  for  the  XDS  Sigma  7  computer  at  M.S.U. 
The  associated  statistical  theory  for  this  type  of  model  can  be  found 
in  [Malinvaud,  Chapter  9]. 

The  Winter  Wheat  Equation 

In  fitting  a  polynomial  yield  surface,  the  shape  of  the  surface 
frequently  does  not  conform  to  a  priori  assumptions.   For  example,  there 
can  easily  be  local  areas  where  the  surface  is  not  concave  downward, 
especially  areas  considerable  distance  from  the  mean  point  of  the  indepen- 
dent variable.   This  problem  was  encountered  in  statistically  fitting  the 
winter  wheat  equation. 

In  particular,  the  surface  was  not  concave  with  respect  to  variation 
on  nitrogen  when  total  water  was  bi' 1 ow  about  7.3  inches.   ThLs  condition 
gave  nonsense  results  for  an  optimal  nitrogen  fertilization  decision  rule 
at  very  low  levels  of  soil  water  at  seeding  time.   This  deficiency  was 
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overcome  by  imposing  a  linear  constraint  on  the  parameter  estimates. 
The  second  derivative  of  (2.1)  with  respect  to  N  is 

(2.9)  — ^  =  2(b,+  3b^N  +  b^W). 
3n2      2     3     8 

The  following  constraint  was  imposed, 

(2.10)  b^  +   5bg  =  0 

Referring  to  (2.9),  it  is  seen  that  when  N  =  0,  the  constraint  forces 

7      9  3 

S  y/3N   =  0  when  W  =  5.   This  constraint  worked  well  when  the  N  -term 
was  deleted  from  the  equation  (b-,  =  0).   Since  the  estimate  of  b   was 
only  marginally  significant,  and  the  positive  estimate  of  h-,   seemed  to 
merely  "flatten  out"  the  response  with  respect  to  nitrogen  at  high  levels 
of  nitrogen,  this  term  was  deleted  in  the  final  model.   Considerable 
experimentation  with  b   included  suggested  that  an  additional  constraint 
besides  (2.10)  would  be  required  if  b_  were  left  in  the  equation;  thus 
b^  =  0  is  a  simple  choice  for  that  second  constraint. 

The  final  estimate  of  the  yield  response  equation  is  given  in 
Table  1.   The  coefficient  of  multiple   determination  is  not  very  meaning- 
ful in  this  situation  because  of  the  correlation  across  treatments; 
however  its  value  is  0.86.   Tte   correlation  across  treatments  is  0.66 
which  is  obviously  not  trivial. 

It  is  easily  shown  that  when  the  two  error  components  in  (2.5)  are 
uncorrelated  (a  plausible  assumption  in  this  application)  ,  the  correlation 
p  is  the  ratio  of  the  year  to  year  variance  to  the  total  variance,  which 
is  0.66  in  our  results  here.   Therefore,  about  two-thirds  of  the  unexplained 
variance  in  the  logarithm  of  yields  is  due  to  the  random  shocks  associated 
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with  different  years.   This  relationship  holds  approximately  in  the 
natural  units  of  bushels  per  acre,  also. 

The  parameter  estimate  for  extra  nitrogen  mineralized  during  the 
summer  fallow  year  is  21  pounds,  which  might  seem  a  little  low  based 
on  work  with  chemical  tests.   However,  the  statistical  precision  leaves 
something  to  be  desired,  note  the  standard  error  of  8.3  pounds.   The 
"true  parameter"  could  easily  be  30  (jr  35  pounds  and  be  consistent 
with  the  statistical  results.   I'his  particular  parameter  estimate  seemed 
a  little  sensitive  to  the  model  specification  too.   One  reassuring 
result  is  that  the  barley  yield  equation  gave  essentially  the  same 
estimate  for  this  parameter. 

Significance  levels  implied  by  the  t-ratios  are  only  approximate 
for  two  reasons:   (1)  the  parameters  enter  the  least  squares  criterion 
in  a  nonlinear  way;  and,  (2)  the  correlation  parameter  which  enters  in 
the  generalized  least  squares  criterion  is  estimated  simultaneously  with 
the  parameters.   A  few  analyses  on  what  appeared  to  be  marginally  signi- 
ficant parameters  at  the  5  percent  levels,  suggested  that  the  parameter 
estimates  are  typically  more  significant  than  indicated  by  the  t-ratios. 

A  likelihood  ratio  test  can  be  performed  by  doing  an  analysis  with 
and  without  a  parameter  (or  set  of  parameters)  included.   This  test 

showed  the  fallow  effect  to  be  marginally  significant  at  the  0.02  level, 

2 
and  the  NW   term  was  marginally  significant  at  the  0.03  level.   Both 

of  these  levels  are  smaller  than  suggested  by  the  t-ratios  reported  in 

Table  1. 
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Table  1 


Statistical  Results  for  Winter  Wheat 
Response  Surface 


Polynomial 

Regression 

Standard 

Term 

Coefficient 

Error 

t-ratio 

constant 

-2.35 

5.22 

-0.451 

N 

-0.417 

0.0713 

-5.84 

N^ 

0.00168 

0.000397 

4.24 

W 

3.72 

0.992 

3.75 

W^ 

-0.140 

0.0427 

-3.28 

NW 

0.0838 

0.0145 

5.77 

N^W 

-0.000336 

0.0000794 

-4.24 

NW^ 

-0.00110 

0.000563 

-1.95 

Fallow  Effect 

6.61 

2.97 

2.23 

Mineralized 

Nitrogen 

on  Fall< 

3w: 

21.1 

8.27 

2.56 

Correlation  among  treatments  =  0.658 
Standard  error  of  estimate  =  0.209 


/, 


(unit  is  natural  logarithm  of  yield—  ) 
46  degrees  of  freedom 


a/ 

—  A  comparable  measure  in  the  original  unit  of  measure,  at  mean  yield  in 

the  sample,  is  approximately  8.0  bushels  per  acre.   (This  is  one-half 

of  a  68%  confidence  interval  in  the  original  units.) 
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The  fallow  effect  is  a  net  effect  in  addition  to  benefits  from 
water  storage  and  extra  nitrogen  mineralized  in  the  soil.   Agronomists 
have  conjectured  that  this  net  effect  is  largely  a  result  of  improved 
stand  when  planting  on  fallow  as  compared  to  seeding  in  stubble. 

A  tabular  presentation  of  the  fitted  yield  response  surface  is  given 
in  the  Appendix.   The  interested  reader  can  discern  the  general  structure 
of  the  surface  from  this  table.   A  profit  function  is  also  given  for 
wheat  price  at  $3  per  bushel  and  nitrogen  price  at  $0.20  per  pound. 

It  is  easy  to  derive  an  optimal  decision  rule  for  nitrogen  fertiliza- 
tion of  wheat.   By  an  optimal  decision  rule,  we  mean  a  schedule  showing 
the  amount  of  nitrogen  fertilizer  to  apply  for  an  observed  level  of  soil 
water  at  seeding  time  such  that  conditional  expected  profit  is  maximized. 
An  equation  with  N  the  dependent  variable  and  W  the  independent  variable 
is  obtained  by  solving  for  N  in  the  following  equation, 

(2.11)  P  -^v  =  p 

W  "dN    N 

where  P^^  and  P   are  prices  of  wheat  and  nitrogen;  respectively,  and 
9y/9N  is  the  slope  of  the  yield  response  surface  with  respect  to  nitrogen 
while  holding  water  fixed.   Solution  of  (2.11)  using  the  notation  of 
(2.1)  gives 

(2.12)  N  =  (Pfj/I^w  -  ^1  -  'V^W  -  l'8W^)/2(b2  +  b7W)  . 

The  nc-xt  step  is  to  make  a  change  of  variable  from  W,  which  is  total 
water,  to  say  S,  wliich  is  soil  water.   This  is  done  by  simply  defining  W 
as  the  sum  of  soil  water  and  mean  precipitation  for  May  through  July. 
Therefore,  we  use  the  substitution 
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(2.13)   W  =  S  +  6.6 

in  (2.12)  since  the  long-term  mean  is  6.6  inches  for  May-July  precipita- 
tion in  the  area  of  the  experiment. 

As  can  be  seen  from  (2.12),  only  the  price  ratio ,  P,,/P  ,  is  important 

N   W 

in  determining  optimal  nitrogen  fertilization.   The  decision  rule  is  given 

in  Table  2  for  $3  wheat  and  $0.20  nitrogen  (P.,/P,,  =  .0667).   It  is  readily 

N   W 

seen  from  the  table  that  the  decision  rule  is  relatively  constant  for 

soil  water  in  e>'cess  of  3  inches.   The  results  do  not  change  much  with  wheat 

prices  except  at  soil  water  levels  below  3  inches. 

An  examination  of  the  profit  function  in  Appendix  Table  1  will  show 
that  profit  is  quite  "flat"  with  respect  to  variations  on  nitrogen  in 
the  vicinity  of  maximum  profits.   Therefore,  an  error  of  less  than  10 
pounds  from  the  optimum  rate  of  nitrogen  application  will  have  a  trivial 
cost  in  lost  profits.   One  might  conclude  that  using  Table  2  without  regard 
to  prices  would  be  a  reasonable  approximation  to  optimal  fertilization. 
Since  profit  variability  tends  to  increase  with  nitrogen  fertilization, 
it  would  be  better  to  error  on  the  low  side  with  respect  to  applied 
nitrogen. 

The  Barley  Yield  Equation 
The  barley  equation  showed  declining  yields  with  respect  to  total 
water  at  very  low  levels  of  water  under  the  initial  estimation.   A  con- 
straint was  introduced  which  forced  incremental  yield  with  respect  to 
water  to  be  zero  at  the  lower  bound  of  the  experimental  data.   The  barley 
response  surface  was  free  from  any  apparent  anomalies  when  this  constraint 
was  applied. 
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Table  2 

a/ 
Nitrogen  Fertilizer  Decision  Rule  for  Winter  Wheat- 


Soil  Water 


Nitrogen  Fertilizer  (lbs/acre) 


(inches  plant  available)  Stubble  Fallow 

1.0  52  32 

1.5  59  39 

2.0  64  44 

2.5  68  48 

3.0  71  51 

3.5  73  53 

4.0  75  55 

4.5  76  56 

5.0  77  57 

5.5  77  57 

6.0  78  58 

6.5  78  58 

7.0  78  58 

7.5  78  58 

8.0  78  58 

8.5  78  58 

9.0  78  58 

9.5  78  58 

10.0  78  58 

10.5  77  57 

11.0  77  57 

11.5  76  56 

12.0  76  56 

12.5  76  56 

13.0  75  55 


Prices  are  $3  wheat  and  $0.20  nitrogen 
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The  empirical  equation  and  statistical  results  are  given  in  Table  3. 
The  low  t-ratios  on  the  interaction  terms  NW ,  N  W,  and  NW^  suggested  a 
model  with  only  NW  as  an  interaction.   Results  of  tliis  equation  are  given 
in  Table  4.   A  comparison  of  results  in  Tables  3  and  4  would  suggest 
thai  I  he  two  inliT.iil  ions  N'^W  .ind  NW^  could  br  dc  1  (.■  trcl  ,  but  a  ILkeliliood 
ratio  tcsL  i  ml  i  iMLcd  Lliat.  Lhesf  Iwo  ic^rins  L<>K<'l-bcr  arf  marginally  signi- 
ficant at  tiie  JO  percent  level.   The  decision  was  made  to  use  the  results 
of  Table  3,  with  all  the  interaction  terms,  in  the  economic  analysis. 

The  standard  error  of  the  estimate  is  much  higher  for  barley  than 
wheat,  indicating  a  poorer  explanation  of  barley  yield  variation.   This 
is  also  suggested  by  the  coefficient  of  multiple  determination  which  is 
0.59  for  barley  compared  to  0.86  for  wheat.   Nevertheless,  the  allocation 
of  unexplained  variance  between  treatments  and  year  to  year  shocks  is 
about  equal  for  the  two  crops  since  the  correlation  across  treatments 
is  0.64  for  barley  and  0.66  for  winter  wheat. 

The  net  effect  of  planting  on  fallow  instead  of  stubble  is  considerably 
greater  for  barley  than  winter  wheat,  14.7  bushels  per  acre  instead  of  6.6. 
Estimated  amount  of  nitrogen  mineralized  during  the  year  of  fallow  is 
nearly  the  same  as  estimated  by  the  winter  wheat  equation,  19  versus  21 
pounds  per  acre,  which  is  reassuring. 

A  decision  rule  for  nitrogen  fertilization  of  barley  is  given  in  Table 
5  for  a  barley  price  of  $1.90  per  bushel  and  nitrogen  price  of  $0.20  per 
pound.   As  in  the  case  of  winter  wheat  fertilization,  the  profit  function 
is  quite  "flat"  with  respect  to  variations  in  nitrogen,  holding  water 
constant.   Therefore,  little  will  be  lost  by  using  the  rule  in  Table  5 
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Table  3 

Statistical  Results  for  Barley  Yield 
Response  Surface 


Polynomial 
Term 


Regression 

Standard 

Coefficient 

Error 

t-rat 

99.0 

19.5 

5.07 

0.208 

0.407 

0.510 

0.00106 

0.00278 

0.381 

-  24.8 

6.43 

-3.87 

constant 

N 

W 


w 


2.47 


0.634 


3.90 


W^ 

-  0,0724 

0.0183 

-3.95 

NW 

-   0.0317 

0.0577 

-0.549 

N  W 

-   0.000314 

0.000238 

-1,32 

NW 

0.00299 

0.00209 

1,43 

Fallow 

Effect 

14.7 

3.45 

4,26 

Minera; 

Lized 

Nitrogen 

on 

Fallow: 

18.8 

10.9 

1,72 

Correlation  among  treatments  =  0.642 
Standard  error  of  estimate  =  0,304   , 
(unit  is  natural  logarithm  of  yield—  ) 

46  degrees  of  freedom 


a/ 

—  A  comparable  measure  in  the  original  unit  of  measure,  at  mean  yield 

in  the  sample,  is  13,1  bushels  per  acre. 
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Table  4 


A]  ItTiiat  i vo  Modul  for  Barley  Yield 
Response  Surface 


Polynomial 

Term 

Regression 
Coefficient 

Standard 
Error 

t-ratio 

constant 

81.9 

13.7 

5.96 

N 

0.120 

0.116 

1.03 

N^ 

-  0.00247 

0.000993 

-2.49 

W 

-19.6 

4.50 

-4.35 

W^ 

1.93 

0.437 

4.43 

w^ 

-  0.0558 

0.0122 

-4.56 

NW 

0.0137 

0.00650 

2.11 

Fallow 

Effect 

14.5 

3.29 

4.41 

Mineralized 
on  Fallc 

Nitrogen 
3w: 

18.5 

13.4 

1.38 

Correlation  among  treatments  =  0.645 
Standard  error  of  estimate  =  0.308 


48  degrees  of  freedom 
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Table  5 


Nitrogen  Fertilizer  Decision  Rule  for  Barley^' 


Nitrogen  Fertilizer  (lbs/acre) 

Soil  Water 
(inches  plant  available)  Stubble  Fallow 


1.0  13  0 

1.5  14  0 

2.0  16  0 

2.5  17  0 

3.0  19  0 

3.5  21  1 

4.0  23  3 

4.5  25  5 

5.0  27  7 

5.5  29  9 

6.0  31  11 

6.5  33  13 

7.0  35  15 

7.5  37  17 

8.0  39  19 

8.5  42  22 

9.0  44  24 

9.5  46  26 

10.0  48  28 

10.5  51  31 

11.0  53  33 

11.5  55  35 


a/ 

—  Prices  are  $1.90  barley  and  $0.20  nitrogen 
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even  when  prices  do  not  correspond;  remember  it  is  the  price  ratio  P  /P 

N   W 

which  determines  the  optimal  rate  of  nitrogen.   A  price  of  $1.93  was  the 
ten-year  average  [irice  lor  barU-y  in  coiislanL  1976  dollars  for  tlir  period 
1966-7S,  and  the  counterpart  lor  winu-r  wheat  was  $3.04  per  husliel. 
These  are  average  market  prices  received  by  farmers  in  Montana. 

Analysis  of  Cropping  Strategies 

First  of  all,  we  must  devise  some  method  to  make  the  experimental 
yield  equations  comensurate  with  experienced  farm  yields.   Proportional 
adjustments  of  the  equations  would  seem  to  make  the  most  sense.   A  recent 
study  of  production  costs  for  Chouteau  County,  based  on  information  from  a 
panel  of  progressive  farmers,  gave  panel  estimates  of  35  and  40  bushels 
per  acre  for  winter  wheat  and  barley,  respectively.   These  yields  are 
used  to  determine  a  factor  for  adjustment.   Both  yields  are  for  crops 
planted  on  fallow. 

Mean  yields  for  planting  on  fallow  are  estimated  from  the  experimental 
yield  equations  by  averaging  predicted  yields  for  each  year  of  the 
experiment.   Predicted  yields  are  calculated  by  determining  the  optimal 
amount  of  nitrogen  fertilizer  to  apply  for  observed  soil  water  at  seeding 
time  during  a  given  year,  and  then  taking  total  water  as  the  sum  of  soil 
water  and  mean  precipitation  for  May  through  July.   This  procedure 
provides  a  longer  period  than  the  experiment  for  averaging  across  seasonal 
precipitation  levels,  but  must  rely  on  only  7  years  of  data  for  plant 
available  soil  water  after  a  year  of  fallow.   Some  improvement  could  be 
made  in  the  procedure  by  fitting  a  frequency  curve  to  soil  water  measure- 
ments, but  it  was  considered  rather  questionable  in  view  of  the  extra 
computation. 
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Average  yields  obtained  by  the  above  procedure  for  the  fitted 
equations  were  53.6  and  61.9  bushels  per  acre  for  winter  wheat  and  fallow, 
respectively.   If  we  adjust  winter  wheat  yield  downward  by  a  factor  so  that 
yield  equals  35  bushels,  the  factor  is  0.65,  i.e.,  (0.65)(53.6)  =  35. 
Rather  surprisingly,  a  comparable  adjustment  on  barley  yield  gives  the  same 
factor  (to  two  significant  digits). 

The  yield  equations  applicable  for  farm  yields  are  as  given  in 
Tables  1  and  3  for  winter  wheat  and  barley,  respectively,  except  that  each 
equation  is  multiplied  by  0.65.   In  the  economic  analysis  to  follow,  we 
can  adjust  price  of  the  product  downward  by  the  factor  0.65  and  use  the 
experimental  yield  equations  directly  because  gross  revenue  is  adjusted 
the  same  either  way. 

Soil  Moisture  Analysis 

We  need  an  estimate  of  plant  available  soil  water  frequencies  at 
seeding  time  in  the  fall.   The  only  data  are  for  seven  years  during  the 
experiment;  the  mean  for  that  period  was  4.23  inches  and  the  standard 
deviation  was  2.59.   A  symmetric  distribution  like  the  normal  is  inappro- 
priate because  of  truncation  at  zero,  less  than  two  standard  deviations 
from  the  mean.   We  also  have  truncation  at  field  capacity,  but  to  five  feet 
for  this  soil;  the  upper  bound  is  about  1  >  inches  which  is  considerably 
farther  from  the  mean  than  the  lower  bound. 

The  small  number  of  observations  hardly  warranted  an  elaborate 
evaluation  of  various  functional  forms;  the  power  of  any  test  for  goodness 
of  fit  would  be  extremely  low.   Essentially  a  log-normal  distribution  was 
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used  with  truncation  at  the  upper  bound  of  soil  water  capacity.   A  small 
change  from  the  truncated  log-normal  distribution  was  made  to  better  deal 
with  the  small  sample  available. 

After  a  logarithmic  transformation  (natural  base) ,  mean  plant 
available  water  was  1.22  and  the  standard  deviaLion  was  0.795.   Let  x  be 
any  plant  available  soil  water  measuri'iiient  drawn  at  random,  and  suppose 
the  mean  of  x  is  |i   and  s   is  the  usual  sample  estimate  of  the  standard 

X         X 

deviation  of  x.   Then,  if  x  is  normally  distributed, 

X  -  u 

X 

=  t 


s 

X 

is  distributed  according  to  the  t-distribution  with  degrees  of  freedom 
one  less  than  sample  size.   For  x  defined  as  the  natural  logarithm  of 
plant  available  soil  water,  we  use 

.n  -,,x     X  -  1-22 
(2-1^)      0.795    ^  ^  ' 

and  assume  that  the  t-distribution  with  six  degrees  of  freedom  is  a  good 
approximation.   If  x  were  normal,  the  sample  mean  x  has  been  used  in  place 
of  the  actual  mean,  p  ,  as  an  approximation.   We  used  (2.14)  instead 
of  the  standardized  normal  because  the  small  number  of  degrees  of  freedom 
gives  a  very  imprecise  estimate  of  the  standard  deviation,  and  it  was  thought 
the  t-distribution  fitted  as  described  would  be  better.   A  tabular 
approximation  to  the  fitted  distribution  is  given  in  Table  6. 

A  comparable  procedure  is  required  to  estimate  the  probability  distri- 
bution for  plant  available  soil  water  in  the  spring  after  a  crop  of  winter 
wheat.   The  mean  of  soil  water  was  5.8  inches,  and  the  standard  deviation 
was  2.19,  for  the  spring  data  at  seeding  time  (same  seven  years).   Notice 
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Table  6 

Probability  Distribution  for  Autumn  Soil  Water 
After  a  Crop  of  Barley 


Water 
Interval 
(inches) 

Midpoint 
(inches) 

Probability 

Cumulative  Probability 
(upper  bound  of  interval) 

0-0.5 

0.25 

.0284 

.0284 

0.5-1.0 

0.75 

.0660 

.0944 

1.0-1.5 

1.25 

.0914 

.1858 

1.5-2.0 

1.75 

.1005 

.2863 

2.0-2.5 

2.25 

.0988 

.3851 

2.5-3.0 

2.75 

.0906 

.4757 

3.0-3.5 

3.25 

.0795 

.5552 

3.5-4.0 

3.75 

.0681 

.6233 

4.0-4.5 

4.25 

.0584 

.6817 

4.5-5.0 

4.75 

.0490 

.7307 

5.0-5.5 

5.25 

.0419 

.7726 

5.5-6.0 

5.75 

.0352 

.8078 

6.0-6.5 

6.25 

.0298 

.8376 

6.5-7.0 

6.75 

.0253 

.8629 

7.0-7.5 

7.25 

.0217 

.8846 

7.5-8.0 

7.75 

.0187 

.9033 

8.0-8.5 

8,25 

.0173 

.9206 

8.5-9.0 

8.75 

.0129 

.9  335 

9.0-9.5 

9.25 

.0123 

.9458 

9.5-10.0 

9.75 

.0109 

.9567 

lO.O-Ll.O 

J0.5 

.0178 

.9745 

11.0-12.0 

11.5 

.0142 

.9887 

12.0-13.0 

12.5 

.0113 

1.0000 
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that  the  mean  is  larger  and  the  standard  deviation  smaller  than  for  the 
fall  measurements.   Measurements  in  the  spring  are  only  to  four  feet 
instead  of  five  feet  as  in  the  fall.   Total  capacity  of  tlie  soil  to 
four  feet  is  about  1.  inches.   Notice  that  the  sample  mean  of  5.8  is 
ai)Oiit  half  way  be  I  ween  Lhc  Lwo  {pounds. 

In  view  of  the  apparent  symmetry  for  soil  water  measurements  in  the 
spring,  we  did  not  transform  the  data  to  logarithms,  but  applied  a  truncated 
t-distribution  to  the  raw  data.   Truncation  was  at  both  ends  of  the  data, 
and  arbitrarily  chosen  sjonmetrically  around  the  mean  with  a  total  range  of 
10  inches.   Thus  the  soil  water  variable  is  confined  to  the  interval 
5.8+5  inches.   A  tabular  approximation  to  the  distribution  is  given  in 
Table  7. 

Basic  Trade-off  Equation 

A  formal  mathematical  optimization  model  can  be  used  to  evaluate 
various  cropping  strategies  (see  Burt  and  Allison,  or  Burt  and  Johnson), 
but  we  use  an  approximation  here,  which  in  this  particular  application, 
turns  out  to  give  the  same  results  as  more  elaborate  methods  would  have. 
One  reason  for  the  simplification  is  that  soil  water  is  essentially 
independent  from  year  to  year  at  the  experimental  site;  that  is  to  say, 
the  soil  water  measurement  this  fall  has  no  influence  on  the  same  measure- 
ment next  fall.   A  similar  result  holds  for  spring  soil  water  measure- 
ments.  Simple  correlations  on  water  measurements  in  adjacent  years  for 
each  spring  and  fall  gave  negative  correlations  which  were  essentially 
zero,  giving  no  sujiport  for  a  [positive  correlation  as  expected.   'I'here  was 
a  positive  correlation  between  spring  soil  water  and  fall  soil  water  of 
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Table  7 


Probability  Distribution  for  Spring  Soil  Water 
After  a  Crop  of  Winter  Wheat 


Water 

Interval 

Midpoint 

Probability 

(inches) 

(inches) 

0.8-1.3 

1.05 

.0123 

1.3-1.8 

1.55 

.0171 

1.8-2.3 

2.05 

.0233 

2.3-2.8 

2.55 

.0313 

2.8-3.3 

3.05 

.0412 

3.3-3.8 

3.55 

.0530 

3.8-4.3 

4.05 

.0654 

4.3-4.8 

4.55 

.0772 

4.8-5.3 

5.05 

.0871 

5.3-5.8 

5.55 

.0921 

5.8-6.3 

6.05 

.0921 

6.3-6.8 

6.55 

.0871 

6.8-7.3 

7.05 

.0772 

7.3-7.8 

7.55 

.0654 

7.8-8.3 

8.05 

.0530 

8.3-8.8 

8.55 

.0412 

8.8-9.3 

9.05 

.0313 

9.3-9.8 

9.55 

.0233 

9.8-10.3 

10.05 

.0171 

10.3-10.8 

10.55 

.0123 

Cumulative  Probability 
(upper  bound  of  interval) 


.0123 
.0294 
.0527 
.0840 
.1252 
.1782 
.2436 
.3208 
.4079 
.5000 
.5921 
.6792 
.7564 
.8218 
.8748 
.9160 
.9473 
.9706 
.9877 
1.000 
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0.58,  but  our  analysis  is  unaffected  by  this  relationship. 

One  side  of  the  trade-off  equation  is  average  annual  net  returns 
from  an  acre  of  land  when  an  optimal  strategy  is  followed.   As  an  approxi- 
mation for  this  value,  we  use  average  annual  net  returns  from  the  alternate 
winter  wheat-fallow  rotation.   If  it  should  turn  out  that  this  rotation 
is  superior  to  any  cropping  strategy  based  on  soil  water  at  seeding  time, 
then  our  approximation  is  perfect. 

Fixed  costs  do  not  enter  into  the  decision  in  making  year  to  year 

evluations  of  which  crop  to  plant,  so  we  deduct  only  variable  costs  to 

get  net  returns.   The  following  notation  is  introduced: 

m  =  plant  available  soil  water  (inches) 

P   =  price  of  winter  wheat  per  bushel  (dollars) 
w 

P   =  price  of  fertilizer  nitrogen  per  pound  (dollars) 

c (m)  =  a  function  describing  the  cost  of  fertilizer  nitrogen  when 

planting  on  stubble  (dollars  per  acre) 
z  =  variable  costs  of  production,  except  nitrogen  fertilizer  costs 

are  excluded  (dollars  per  acre) 
C   =  variable  costs  of  a  year  in  summer  fallow  (dollars  per  acre) 
y(m)  =  conditional  expected  yield  of  winter  wheat  planted  on 
stubble  when  plant  available  soil  water  is  m  inches  at 
seeding  (bushels  per  acre) 
y   =  average  yield  of  winter  wheat  when  planted  on  fallow  (bushels 

per  acre) 
c,  =  average  cost  of  nitrogen  fertilizer  when  planting  on 
fallow  (dollars  per  acre) 
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Nitrogen  costs,  c(m)  and  c,,  are  functions  of  the  price  ratio  P  /P   since 

I  n  w 

optimum  nitrogen  fertilization  is  assumed,  but  for  simplicity  in  the  notation, 
the  price  ratio  was  not  included  in  the  functional  forms.   Also,  the  price 
ratio  is  implicit  in  yields  y(m)  and  y  . 

Elementary  reasoning  would  tell  us  that  when  soil  water  is  high 
enough  in  the  fall  that  expected  net  returns  from  planting  a  crop  of 
wheat  exceeds  average  annual  net  returns  from  the  best  strategy,  then 
we  should  plant  the  wheat  and  vice  versa.   The  condition  where  we  would 
be  indifferent  is  when  the  expected  net  returns  from  planting  are  just 
equal  to  average  net  returns  under  the  best  strategy.   We  are  using  alter- 
nate winter  wheat-fallow  as  an  approximation  for  the  best  strategy. 

Therefore,  our  trade-off  equation  is 

(2.15)  P^^y(m)  -  c(m)  -  z  =  [P^7^  _  z  -  C^  -  I^]/2 

The  left  hand  side  of  (2. 15)  is  expected  net  returns  from  planting  a  crop 
of  winter  wheat  when  observed  soil  water  is  m.   The  right  hand  side  is 
average  annual  net  returns  in  alternate  winter  wheat-fallow;  division  by 
two  allows  for  a  crop  every  other  year. 

We  transpose  z  on  the  left  to  the  right  side  and  simplify  to  get 

(2.16)  P^y(m)  -  c(m)  =  z/2  +  [P^7^  '   ^f   '   ^f]/2 

The  above  equation,  with  m  as  the  unknown  variable,  is  solved  to  determine 

the  threshold  level  of  soil  water  above  which  winter  wheat  is  planted  on 

stubble,  and  below  which,  a  year  of  summer  fallow  is  decided  upon.   This 

result  is  specific  to  given  prices,  P   and  P  . 

w      n 

The  above  discussion  would  lead  one  to  assume  that  we  are  considering 
planting  winter  wheat  on  winter  wheat  stubble,  but  this  is  not  the  case.   The 
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previous  crop  would  have  to  have  been  barley  or  some  other  crop  like 
safflower  that  left  a  light  enough  stubble  to  make  planting  winter  wheat 
feasible  in  this  climate. 

The  analysis  was  done  with  the  following  estimates  of  costs  and  price 
of  nitrogen: 

z  =  $20 

Cf  =  $5 

P   =  $0.20 
n 

The  remaining  terms  in  (2. 16)  depend  on  the  winter  wheat  price  specified. 
Variable  costs  of  production  were  taken  from  (Schaefer,  etal.  1976c)  adjusted 
by  deleting  the  following  items  from  variable  costs:   (1)  nitrogen  fertilizer; 
(2)  crop  insurance;  and,  (3)  labor.   Variable  costs  of  fallow  were  taken 
from  the  same  source  with  labor  costs  deducted. 

The  analysis  for  a  spring  decision  on  whether  to  plant  barley  after  a 
year  with  the  land  in  crop  is  essentially  the  same.   The  only  difference 
is  in  the  interpretation  of  the  left  hand  side  of  (2.16);  y(m)  is  then 
barley  yield  and  P   is  replaced  with  barley  price.*   It  is  assumed  that  costs 
of  production,  z,  are  the  same  for  barley  as  wheat.   Also,  the  nitrogen 
fertilizer  cost  function,  c(m),  must  be  reinterpreted  as  that  appropriate 
for  optimum  barley  fertilization. 

"Optimum"  nitrogen  fertilization  was  determined  with  the  scaled  down 
yield  equation  which  makes  the  rates  lower  than  given  in  Tables  2  and  5  for 
winter  wheat  and  barley,  respectively.   Also,  the  analysis  was  done  for 
several  prices  of  wheat  and  barley. 


*See  page  38-a. 
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*After  the  following  analysis  was  completed,  it  was  realized  that  the 
conditional  expected  net  returns,  given  the  land  is  in  barley  stubble  in 
the  fall,  is  greater  than  average  annual  returns  from  alternate  wheat- 
fallow  by  a  significant  amount  (about  $5  per  acre  with  wheat  at  $3  per 
bushel).   In  view  of  this  fact,  a  more  precise  trade-off  equation  would 
be  to  change  the  left  hand  side  of  (2.15)  to  the  average  of  net  returns 
from  barley  (given  soil  moisture)  and  expected  net  returns  from  the 
following  year  if  the  land  in  the  autumn  were  in  barley  stubble.   This 
two-year  average  would  introduce  the  extra  advantage  which  planting 
barley  has  on  expected  net  returns  as  viewed  in  the  fall  after  the  barley 
has  been  harvested.   The  difference  in  results  is  quite  trivial  as  far 
as  estimating  an  optimal  decision  rule  based  on  an  expected  profit  criterion, 
but  the  measure  of  expected  opportunity  loss  used  below  to  estimate  costs 
of  deviating  from  the  "otpimal"  strategy  is  significantly  affected.   This 
difference  is  quantified  at  the  end  of  the  report. 
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Relatlve  Price  of  Barley  to  Wheat 

Economic  strategies  for  cropping  decisions  obviously  depend  critically 
on  relative  prices  of  barley  and  wheat.   A  regression  was  run  with  barley 
price  the  dependent  variable  and  winter  wheat  price  the  independent 
variable  (an  arbitrary  choice).   In  the  raw  data  for  ten  years  (1966-75), 
using  average  prices  received  by  Montana  producers,  the  correlation  was 
0.978  which  suggests  a  strong  systematic  relationship  between  the  two 
prices.   However,  general  changes  in  the  price  level  tend  to  exaggerate 
the  correlation  in  prices. 

Consequently,  the  prices  were  deflated  by  prices  paid  by  farmers 
for  production  items,  and  the  correlation  was  0.958,  still  quite  substan- 
tial.  Measured  in  1976  dollars,  the  standard  error  of  the  estimate  for 
barley  price,  given  winter  wheat  price,  was  $0.18.   The  fitted  regression 
equation  was 

(2.17)   Pjj  =  0.455  +  0.484?^   , 

where  Pi^  and  P^^  are  barley  and  winter  wheat  prices,  respectively.  Mean 
prices  (1976  dollars)  were  $3.04  and  $1.93  for  winter  wheat  and  barley, 
respectively.   Selected  pairs  of  prices  from  (2.17)  are  given  below. 

Table  8 
Relative  Prices  of  Wheat  and  Barley 


Pw 

2.50 

3.00 

3.50 

4.00 

4.50 

5.00 

h 

1.67 

1.91 

2.15 

2.39 

2.63 

2.88 
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Relative  prices  may  deviate  considerably  from  the  relationship  of 
(2.17),  but  for  planning  purposes  and  economic  analysis  such  as  being 
pursued  here,  prices  in  the  above  table  should  be  a  good  guide.   The 
analysis  is  performed  for  the  above  prices  from  $2.50  to  $4.00  for  winter 
wheat. 

Criterion  of  Maximum  Average  Profits 

The  trade-off  equation  (2.16)  was  solved  for  each  of  the  four  wheat 
prices  and  the  results  are  summarized  in  Table  9.   These  results  constitute 
the  "optimal"  decision  rule  for  decisions  in  the  fall  when  the  land  is 
in  barley  stubble  so  that  planting  winter  wheat  is  feasible.   By  "optimal," 
we  mean  within  the  context  of  the  model  and  the  empricial  yield  response 
functions  which  are  being  used  for  the  analysis.   In  actual  farm  circum- 
stances, yields  from  alternate  winter  wheat-fallow  might  be  lower  than  if 
barley  were  worked  into  the  rotation.   These  complications  are  introduced 
later. 

For  the  spring  decision  and  a  criterion  of  maximum  average  profits, 
solution  of  (2.16),  adjusted  on  the  left  hand  side  for  barley,  gives  a 
critical  soil  moisture  level  of  8.9  inches  when  wheat  price  is  $2.50  and 
barley  price  is  $1.67.   At  the  other  pairs  of  prices  in  Table  8  (up  to 
$4  for  wheat  price),  barley  would  never  be  planted.   That  is  to  say,  barley 
would  not  be  superior  to  alternate  winter  wheat-fallow  at  any  observed  soil 
water  level  in  the  spring. 
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Table  9 

Autumn  IJcfision  Rule  Maximizing  Mean  Profits 


Wheat 

Critical  Soi 

Price 

Water 

Leveli. 

(dollars) 

(inches) 

2.50 

4.0 

3.00 

3.5 

3.50 

3.2 

4.00 

3.0 

Right  Hand 
Side  of  Equation 


(dollars) 
46.00 
54.50 
63.25 
71.75 


a/ 


The  decision  rule  is  plant  winter  wheat  if  plant  available  soil  water 
exceeds  the  critical  level  and  summer  fallow  if  it  is  below  the  critical 
level. 
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Expected  Opportunity  Loss  Analysis 

Winter  wheat  yields  can  be  expected  to  increase  when  barley  is 
introduced  into  the  cropping  sequence  because  of  benefits  from  disease 
and  weed  control,  as  well  as  other  more  diverse  agronomic  benefits.   This 
is  particularly  true  under  actual  farm  conditions  as  contrasted  to  the 
more  intensive  controls  imposed  under  experimental  conditions.   Consequently, 
the  rather  simple  comparisons  made  with  the  economic  trade-off  equation  of 
(2.16)  could  easily  result  in  barley  not  being  planted  frequently  enough. 

The  following  analysis  is  directed  at  estimating  the  average  oppor- 
tunity loss  associated  with  various  decision  rules,  where  opportunity  loss 
is  the  average  difference  in  net  returns  under  the  "optimal"  decision  rule 
according  to  (2.16)  and  a  specific  decision  rule  which  plants  barley  with 
a  greater  frequency  and  also  plants  winter  wheat  on  barley  stubble  with  a 
greater  frequency.   The  average  opportunity  loss  can  then  be  evaluated  in 
the  light  of  prospective  winter  wheat  yield  increases  and  control  of  saline 
seep  problems. 

A  decision  rule  is  specified  by  two  numbers:   (1)  the  critical  value 
of  autumn  soil  water  and  (2)  the  critical  value  of  spring  soil  water.   Let 
us  call  these  two  critical  values  M^  and  H„    for  the  autumn  and  spring 
decisions,  respectively.   We  use  the  convention  M  =  co  to  mean  that  barley 
is  never  planted  in  the  spring.   Recall  that  a  crop  is  planted  if  plant 
available  soil  water  exceeds  the  critical  value  and  vice  versa. 

For  any  pair  of  critical  values  (M^,  M^),  we  can  use  the  probability 
Tables  6  and  7  to  obtain  the  probabilities  of  pJanting  crops  in  the  autumn 
and  spring.   For  example,  take  M^  =  2.0  and  M  =  5.3;  from  Table  6,  it  is 
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seen  that  tiie  probabil  Lty  of  nut  |)lanting  winter  wheat  on  stubble  is 
0.2865  (last  column,  line  A);  likewise,  the  probability  of  not  planting 
barley  on  stubble  in  the  spriny  is  O.A079  (last  column  of  Table  7, 
line  9).   The  probabilities  of  planting  crops  are  obtained  by  subtracting 
these  probabilities  from  one. 

The  restriction  is  made  in  the  analysis  here  that  winter  wheat  will 
always  be  planted  after  a  year  of  summer  fallow.   This  specification  is 
warranted  by  a  simple  comparison  of  average  net  returns  from  each  barley 
and  winter  wheat  planted  on  fallow.   An  examination  of  the  paired  prices 
of  winter  wheat  and  barley  in  Table  8  suggests  that  the  higher  prices, 
the  more  advantage  winter  wheat  will  have.   A  comparison  of  mean  net  re- 
turns at  prices  corresponding  to  average  prices  (1976  dollars)  over  the 
period  1966-75  shows  a  difference  of  about  $25  per  acre  in  favor  of  winter 
wheat.   This  result  is  with  yields  scaled  down  from  experimental  yields 
to  35  and  40  bushels  per  acre  for  winter  wheat  and  barley,  respectively. 
Even  at  $2.50  wheat  and  $1.67  barley  prices,  the  difference  is  about  $17 
per  acre.   Of  course,  something  like  malting  barley  prices  could  change 
this  substantially  for  any  given  year  and  producer.   Nevertheless,  this 
evidence  is  rather  convincing  as  a  general  rule. 

The  cropping  systems  strategy  can  be  summarized  as  follows: 

1.  Plant  winter  wheat  on  summer  fallowed  land; 

2.  Plant  winter  wheat  on  barley  stubble  if  plant  available  soil 

water  exceeds  M  :  and, 
a' 

3.  Plant  barley  on  winter  wheat  stubble  if  plant  available  soil 
water  exceeds  Mg. 
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The  state  of  land  use  in  any  particular  year  is  either  fallow,  winter  wheat, 

or  barley. 

Continuing  with  the  example  where  M  =2.0  and  M  =  5.3,  a  set  of 

conditional  probabilities  can  be  defined  for  land  use  next  year  for  a  given 

land  use  this  year.   These  probabilities  are  given  for  this  example  in 

Table  10.   The  first  row  is  obvious  since  winter  wheat  is  always  planted 

on  fallow.   The  probabilities  in  the  second  row  emanate  from  the  probability 

of  soil  water  exceeding  or  not  the  value  M  =5.3.   As  indicated  earlier 

s 

0.4079  is  the  probability  of  soil  water  not  exceeding  5.3  in  the  spring 

(read  from  Table  7) ,  and  is  thus  the  probability  of  a  year  of  fallow  after 

a  year  of  winter  wheat.   Its  complement,  0.5921,  is  the  probability  of 

planting  barley  after  winter  wheat,  and  there  is  no  chance  of  a  second  year 

of  winter  wheat.   The  probability  of  a  year  of  fallow  after  barley  is  the 

probability  of  soil  water  in  the  fall  not  exceeding  M  =2.0  (read  from 

a 

Table  6)  which  was  obtained  earlier  as  0.2863.   Its  complement  is  the  pro- 
bability of  planting  winter  wheat  after  a  year  of  barley,  and  again,  the 
chance  of  a  second  year  of  the  same  crop  is  zero. 

Table  10 

Conditional  Probabilities  for  Land  Use 

With  M  =  2. 0  and  M  =  5. 3 
a  s 

Current 

Land  Use  Probabilities  of  Land  Use  Next  Year 

Fallow  Wheat  Barley 

Fallow  0  1.0  0 


Wheat  0.4079  0  0.5921 

Barley  0.2863  0.7137  0 


Wheat 

1. 

0 

0 

0. 

7137 
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When  land  use  from  year  to  year  obeys  a  probabilistic  pattern  de- 
scribed by  a  table  of  conditional  probabilities,  like  the  example  problem 
in  Table  10,  it  is  called  a  Markov  chain.   By  using  the  elementary  prin- 
ciples of  probability  theory,  the  probability  of  land  use  at  any  future 
date  can  be  calculated  for  any  given  land  use  at  the  present  time.   An 
important  property  of  such  Markov  chains  is  the  existence  of  a  set  of 
equilibrium  probabilities  which  have  the  following  interpretation.   If 
a  year  far  into  the  future  were  picked  at  random,  probabilities  of  land 
use  would  be  independent  of  current  land  use,  and  these  probabilities  are 
unique  as  well  as  easily  calculated  from  a  table  such  as  Table  10  (Kemeny, 
Snell,  and  Thompson). 

These  equilibrium  probabilities  for  a  specified  decision  rule  will 
be  referred  to  as  relative  frequencies  of  land  use  because  they  measure 
the  long-run  relative  frequencies  of  land  use,  regardless  of  current  land 
use.   These  relative  frequencies  are  the  statistical  counterpart  of  the 
proportions  of  time  that  land  is  in  various  uses  under  a  conventional 
rotation.   If  you  please,  think  of  relative  frequencies  of  land  use  as 
describing  a  probabilistic  rotation.   One  word  of  caution,  however,  the 
same  relative  frequencies  could  be  associated  with  more  than  one  cropping 
strategy  (decision  rule)  and  the  economic  consequences  of  these  strategies 
would  not  be  the  same.   In  other  words,  relative  frequencies  of  land  use 
do  not  reveal  all  of  the  information  and  consequences  of  a  particular 
strategy. 

In  the  analysis  of  Markov  chains  for  this  research,  the  relative 
frequencies  were  calculated  by  taking  powers  of  the  matrices  of  probabilities 
(Kemeny,  et  al) .   The  matrices  are  only  three  dimensional  and  n  iterations 
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give  a  matrix  to  the  power  of  2*^,  which  makes  the  method  fast  and  simple 
for  computer  use.   Fewer  tlian  ten  iterations  were  required  to  get  conver- 
gence which  is  accomplished  when  all  the  rows  of  the  computed  matrix  are 
identical. 

Relative  frequencies  for  the  example  given  in  Table  10  are  27 
percent  fallow,  46  percent  winter  wheat,  and  27  percent  barley.   Approxi- 
mately the  same  frequencies  would  prevail  under  a  fixed  four  year  rota- 
tion of  fallow,  winter  wheat,  barley,  winter  wheat,  but  the  implications 
for  soil  water  and  saline  seep  control  would  be  much  different.   The 
decision  rule  associated  with  Table  10  implies  that  fallow  will  be  used 
only  during  the  driest  years  as  reflected  by   soil  water  in  the  fall  and 
spring,  but  the  fixed  four-year  rotation  completely  ignores  soil  water 
at  seeding  time  as  a  criterion  for  fallowing.   Average  deep  water  perco- 
lation through  the  root  zone  would  be  higher  under  the  fixed  rotation 
than  under  the  decision  rule  which  uses  soil  water  as  a  criterion  for 
fallow. 

The  decision  rules  for  the  spring  decision  were  M  =  "°,  except  at 
a  price  of  $2.50  for  wheat  and  $1.67  for  barley.   Even  then,  M  =  8.9 
which  implies  fallow  after  winter  wheat  about  92  percent  of  the  time 
(see  Table  7  last  column  and  fifth  row  up  from  the  bottom).   Recall 
M  =  <«  implies  to  always  tallow  after  winter  wheat.   These  results 
essentially  lock  the  producer  into  winter  wheat-fallow  with  barley  planted 
only  on  rare  occasions  when  prices  are  low  and  soil  water  is  very  high  in 
the  spring. 


-47- 


Conditional  expected  opportunity  losses  can  be  calculated  for 

violating  the  maximum  expected  profit  decision  rule.   We  are  only  interested 

in  violations  which  mean  planting  a  crop  with  "too  little"  soil  water.   For 

example,  let  plant  available  soil  water  be  3.75  inches  in  the  fall  after 

a  crop  of  barley,  and  suppose  wheat  price  is  $2.50.   From  Table  9,  it  is 

seen  that  the  "optimum"  value  of  M   is  A.O,  so  an  observed  level  of  water 

a 

equal  to  3. 75  would  imply  fallowing  if  we  followed  the  decision  rule.   What 

would  be  the  average  opportunity  cost  of  planting  winter  wheat  instead  of 

fallowing? 

Returning  to  our  trade-off  equation  (2.16),  setting  m  =  3.75  will 

result  in  the  left  hand  side  being  smaller  than  the  right  hand  side;  this 

difference  is  conditional  expected  opportunity  loss.   The  computed  difference 

is  $1.44.   Turning  to  Table  7,  note  that  3.75  is  the  midpoint  of  the  soil 

water  interval  3.5  to  4.0  inches  which  occurs  with  probability  0.0681. 

Therefore,  reduction  of  M   from  4.0  to  3.5  would  result  in  an  expected 
'  a 

opportunity  loss  of 

(0.0681)(1.44)  =  $0. 10. 
Observe  the  distinction  between  conditional  and  unconditional  expected 
opportunity  loss;  we  dropped  the  redundant  adjective  unconditional  though. 

By  using  the  procedure  described  above  for  computing  conditional 
expected  opportunity  loss  and  expected  opportunity  loss.  Table  11  was 
constructed  which  provides  a  running  account  of  expected  losses  as  M   is 
systematically  reduced  from  the  values  of  Table  9.   Notice  how  conditional 
expected  opportunity  losses  increase  more  rapidly  than  cumulative  expected 
losses.   In  this  context,  conditional  expected  loss  is  like  an  incremental 
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loss  in  going  to  a  drier  and  drier  decision  rule,  while  cumulative  ex- 
pected losses  are  weighted  average  losses  over  the  interval  of  soil 
water  where  wheat  is  being  planted  in  violation  of  the  decision  rule  in 

Table  9.   If  M  is  chosen  as  the  lower  bound  of  any  interval  in  the  first 
a 

column  of  Table  11,  cumulative  expected  loss  for  that  row  is  expected 
opportunity  loss  associated  with  the  autumn  decision  rule,  and  the  last 
column  of  the  table  is  the  probability  of  planting  wheat  on  stubble. 

As  an  example,  let  us  examine  Table  11  under  a  price  of  $3.00  for 
wheat.   In  the  soil  water  interval  3.5  to  4.0  inches,  there  are  no  con- 
ditional expected  opportunity  losses  because  the  critical  soil  water  level 
under  an  average  maximum  profit  decision  rule  is  3.5  inches  which  lies 
at  the  lower  boundary  of  the  interval;  losses  do  not  occur  unless  a  decision 
rule  is  used  where  wheat  is  planted  with  less  than  3.5  inches  of  plant 
available  water.   In  the  interval  3.0  to  3.5  inches,  conditional  expected 
losses  are  seen  to  be  $2.14  which  is  the  difference  between  $54.50  and 
conditional  expected  net  returns  from  planting  barley  with  plant  available 
water  equal  to  3.25  inches  (midpoint  of  interval).   The  $54.50  used  as  a 
basis  to  calculate  the  loss  is  the  right  hand  side  of  the  trade-off  equation 
(2.16),  which  is  average  net  returns  in  alternate   winter  wheat- fallow. 
The  ]ast  two  rows  for  the  conditional  expected  losses  column  under  $3  wheat 
is  calculated  by  the  same  procedure. 

Notice  how  rapidly  these  conditional  expected  losses  increase  as  the 
decision  rule  is  adjusted  to  lower  and  lower  critical  soil  water  levels, 
reaching  $9.20  for  the  interval  2.0  to  2.5  inches.   This  last  loss  estimate 
means  that  if  winter  wheat  were  planted  on  stubble  when  there  were  only 
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2.25  inches  (midpoint)  of  plant  available  water,  the  average  sacrifice 

in  net  income  per  acre  would  be  $9.20.   Of  course,  this  estimate  imputes 

nothing  to  improved  yields  from  an  integrated  strategy  that  gets  barley 

into  the  rotation  (the  stubble  has  to  be  barley),  nor  is  any  credit 

being  given  to  improved  saline  seep  control;  these  aspects  are  dealt  with 

shortly. 

The  cumulative  expected  loss  column  in  Table  11  is  simply  a  running 

total  of  the  sum  of  the  cross  products  of  the  probability  column  and 

conditional  expected  loss  column.   This  cumulative  expected  loss  column 

can  best  be  interpreted  by  an  example.   Again  using  $3  wheat,  consider  the 

row  for  interval  2.5  to  3.0  inches  of  soil  water  where  cumulative  expected 

loss  is  $0.68  in  Table  11.   If  the  critical  soil  water  level  M  were 

a 

taken  as  2.5  inches,  then  any  year  in  which  barley  stubble  was  available 
for  recropping  would  have  an  average  net  return  of  $0.68  per  acre  below 
what  it  would  have  if  the  "optimum"  critical  value  of  3.5  inches  were 
used.   If  the  critical  soil  water  level  were  pushed  on  down  to  2.0 
inches,  the  average  sacrifice  in  such  a  year  after  barley  would  be  $1.59 
per  acre.   These  expected  opportunity  losses  in  the  cumulative  expected 
loss  column  of  Table  11  are  conditional  on  the  land  being  in  a  stubble  sucli 
as  barley  which  makes  recropping  with  winter  wlieat  feasible.   This 
"conditional"  aspect  must  not  be  confused  with  the  conditional  expected 
loss  columns  of  Table  11  which  imply  a  particular  soil  water  interval  as 
well  as  the  existence  of  barley  or  other  feasible  stubble. 

The  same  kind  of  economic  analysis  is  applicable  to  the  spring  decision 
when  a  decision  must  be  made  on  whether  to  plant  barley  after  winter  wheat 
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or  to  summer  fallow  the  land.   The  counterpart  of  Table  11  is  Table  12 

for  the  spring  decision.   Substantially  more  rows  are  required  because 

the  critical  soil  water  level  under  an  "optimum"  decision  rule  is  relatively 

high;  in  fact,  it  is  field  capacity  for  all  prices  except  $2.50  wheat  where 

it  is  8.9  inches. 

By  using  the  Markov  chain  framework  discussed  earlier  and  Tables  11 

and  12,  we  can  evaluate  various  strategies  for  cropping  decisions.   For 

each  pair  of  critical  soil  water  levels,  M   in  the  fall  and  M^  in  the  spring, 

a  strategy  is  defined.   A  set  of  long-run  relative  frequencies  for  land 

use  exists  for  each  strategy,  and  expected  opportunity  losses  are  available 

in  Tables  11  and  12  for  land  in  barley  stubble  and  winter  wheat  stubble, 

respectively.   Results  are  summarized  in  Table  13. 

Optimal  combinations  of  M  and  M  exist  in  theory  in  the  sense  that 
^  as 

for  a  given  frequency  of  fallow,  expected  opportunity  losses  are  a  minimum. 
These  optimal  combinations  of  M  and  M  would  correspond  to  equal  values 
of  conditional  expected  losses  in  Tables  11  and  12  at  given  prices  of 
grain.   For  example  at  $3  wheat,  M^  =  2.0,  and  Mg  -  6.3;  this  condition 
is  approximately  met  since  conditional  expected  opportunity  losses  are 
$9.20  and  $9.02  in  Tables  11  and  12,  respectively.   From  Table  13,  it  is 
seen  that  this  strategy  gives  a  fallow  frequency  of  33.5  percent  and 
expected  loss  of  $1.24.   We  can  conclude  that  for  this  frequency  of 
fallow  and  $3  wheat,  no  other  strategy  gives  a  smaller  expected  opportunity 
loss. 

However,  serious  limitations  are  encountered  i.-;  applying  this 
principle.   First,  the  results  are  quite  sensitive  to  prices  of  grain,  and 
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second,  the  criterion  often  implies  extremely  small  or  negative  values  for 
M^  at  low  fallow  frequencies  which  are  of  interest.   A  practical  difficulty 
is  considerable  interpolation  in  Tables  11  and  12.   For  these  reasons,  we 
do  not  impose  this  efficiency  criterion  on  the  strategies  examined,  but 
merely  list  several  strategies  in  decreasing  order  according  to  frequency 
of  fallow  (see  Table  13).   Choosing  >I  less  than  2.0  inches  was  considered 

a. 

infeasible  on  practical  grounds  because  of  the  hazards  of  getting  a  stand 
of  winter  wheat.   If  this  constraint  is  accepted  as  valid,  there  is  no 
opportunity  for  application  of  the  efficiency  criterion  when  frequency  of 
fallow  is  less  than  about  33  percent. 

Agricultural  research  indicates  considerable  advantage  in  avoiding 
monoculture  because  of  various  disease  and  weed  problems.  The  practice 
of  alternative  winter  wheat-fallow  has  met  with  considerable  success  in 
spite  of  the  hazards  of  disease  and  weeds,  but  there  also  appears  to  be 
considerable  opportunity  for  increased  yields  by  getting  a  spring  planted 
crop  into  the  cropping  system. 

Bringing  barley  into  the  rotation  would  help  control  cheatgrass  which 
has  a  life  cycle  corresponding  closely  with  winter  wheat.   There  would  also 
appear  to  be  some  advantage  for  the  control  of  wild  oats.   An  occasional 
crop  of  barley  would  tend  to  control  some  of  the  diseases  common  to  winter 
wheat,  especially  soil  borne  diseases. 

Let  us  suppose  for  a  moment  that  bringing  barley  into  the  cropping 

system  with  an  average  frequency  of  20  percent  or  more  would  increase 

average  winter  wheat  yields  by  two  bushel  per  acre,  a  rather  modest  amount. 

In  examining  the  results  of  the  strategy  M  =2.0  and  M   =  5.0  in  Table 

a  s 
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13,  notice  Chat  expected  opportunity  losses  across  all  prices  are  just  a 

little  less  than  the  price  of  a  busnel  of  wheat.   Also,  observe  that  winter 

wheat  is  planted  with  a  relative  frequency  of  46  percent.   Consequently, 

the  two  bushel  increase  in  yield  is  experienced  on  only  46  percent  of  the 

cropland,  on  the  average.   With  $3  wheat,  that  increases  is  equal  to  $2.76 

per  acre  annually  (2  x  3  x  0. 46  =  2.76),  and  the  expected  opportunity  loss 

is  $2.79  (see  Table  13).   Thus  it  is  seen  that  a  two  bushel  increase  in 

average  winter  wheat  yield  attributable  to  barley  in  the  cropping  system 

can  compensate  for  the  expected  opportunity  loss  associated  with  the  strategy 

M^  =  2.0  and  M  =  5.0. 
a  s 

Actually,  the  expected  opportunity  loss  is  somewhat  lower  than  the 
methodology  just  presented  would  lead  us  to  believe.   There  is  a  rather 
subtle  dynamic  element  involved  in  the  trade-off  between  fallowing  and 
planting  barley  in  the  spring  which  is  not  accounted  for  by  the  method 
used  to  calculate  expected  opportunities  loss.   The  additional  complexity 
was  discovered  as  this  final  report  was  being  typed;  time  and  resource 
limitations  did  not  permit  an  entire  new  analysis.   However,  some  limited 
calculations  for  the  $3  wheat  case  showed  expected  opportunity  losses  to 
be  about  one  dollar  for  the  rule  M^  =  2.0  and  M  =5.0.   For  this  strategy, 
it  would  seem  safe  to  claim  that  expected  opportunity  losses  are  probably 
covered  by  a  one  bushel  per  acre  increase  in  winter  wheat  yields  which 
can  be  attributed  to  getting  barley  into  the  rotation  (rather  than  a  two 
bushel  increase  as  suggested  earlier). 

The  measure  of  expected  opportunity  loss  used  in  Table  13  would  be 
appropriate  in  a  somewhat  drier  location  where  land  in  barley  stubble  in 
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the  fall  does  not  have  any  higher  conditional  expected  net  return  than 

average  annual  returns  from  alternate  wheat-fallow.   Consequently,  the 

two-bushel  per  acre  yield  increase  on  winter. wheat  required  to  compensate 

for  expected  opportunity  loss  associated  with  using  the  strategy  M  =2.0 

a 

and  Mg  =  5.0  would  be  appropriate  for  somewhat  drier  winter  wheat  areas 

of  North  Central  Montana.   As  an  informed  guess,  we  might  say  that  for 

areas  in  this  general  region  which  have  14  or  more  inches  of  average  annual 

precipitation,  a  two-bushel  increase  in  winter  wheat  yields  attributable 

to  barley  in  the  cropping  sequence  would  economically  justify  the  strategy 

M_  =  2.0,  M  =  5.0. 
a  s 

Clearly,  the  strategy  M^  =  2.0  and  M  =  5.0  is  quite  "dry"  in  that 
land  is  fallowed  an  average  of  25  percent  of  the  time.   Not  only  is  the 
relative  frequency  of  fallow  low,  but  land  is  not  fallowed  after  barley 
when  soil  water  is  greater  than  2.0  inches  in  the  fall,  and  land  is  not 
fallowed  after  winter  wheat  when  soil  water  in  the  spring  is  greater 
than  5  inches.   It  would  appear  that  such  a  strategy  could  control 
saline  seep  problems  once  the  soil  profile  below  the  root  zone  were 
dried  out  by  a  crop  like  alfalfa,  or  even  safflower  to  some  extent. 

In  cases  where  the  farmer  is  unwilling  to  dry  out  the  deep  subsoil, 
this  strategy  would  probably  control  percolation  be^'ond  the  root  zone,  and 
ultimately  control  the  saline  seep  problem.   But  the  process  would  be  slow, 
just  as  the  problem  evolved  slowly,  and  the  exacerbated  damage  to  the  soil 
around  the  seeps  might  be  an  important  economic  consideration. 

Using  a  fixed  rotation  for  comparison,  fallowing  without  regard  to 
soil  water  will  result  in  summer  fallow  being  carried  out  with  initial  soil 
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water  in  the  fall  greater  than  two  inches  about  70  percent  of  the  time, 
and  in  the  spring,  fallowing  will  take  place  with  soil  water  greater  than 
five  inches  about  two-thirds  of  the  time.   Obviously,  the  saline  seep 
hazards  from  deep  percolation  of  water  are  much  greater  when  soil  water  is 
not  used  as  a  decision  criterion  since  fallowing  is  never  carried  out 
when  soil  water  is  greater  than  two  and  five  inches  in  the  fall  and  spring, 
respectively,  under  the  strategy  M  =  2.0  and  M„  =  5.0. 

The  best  fixed  rotation  for  saline  seep  control  and  relatively  high 
returns  is  four  years  in  length:   fallow-wheat-barley-wheat-   Fallow  only 
occurs  25  percent  of  the  time  and  expected  annual  net  returns  under  $3 
wheat  are  $31.17  compared  to  $34.50  for  alternate  crop-fallow.   This  is  with 
all  variable  costs  deducted,  but  not  fixed  costs.   But  the  strategy  M^  =  2.0, 
M  =5.0  has  an  average  annual  net  return  of  $33.43  and  also  has  a  relative 
frequency  of  fallow  equal  to  25  percent.   This  strategy  has  an  edge  of  $2.25 
per  acre  annually  as  well  as  being  more  effective  in  saline  seep  control. 


III.   THE  SPRING  WHEAT  REGION 

The  research  reported  in  this  section  has  application  in  the  spring 
wheat  areas  of  Northern  Montana  and  North  Dakota.   The  most  important 
empirical  information  required  for  economic  decisions  on  cropping  strategies 
is  a  spring  wheat  yield  response  surface,  where  the  explanatory  variables 
are  soil  water,  precipitation,  and  daily  maximum  temperatures.   Such  an 
equation,  with  sufficiently  broad  application,  can  be  used  as  the  basis 
for  estimating  yields  on  both  sides  of  the  trade-off  equation  (2.16)  in 
Section  II  of  this  report. 

The  decision  problems  are  much  simpler  in  the  spring  wheat  region 
compared  to  a  winter  wheat  area  such  as  was  analyzed  in  Section  II  because 
significant  decisions  on  cropping  are  made  only  in  the  spring.   With  a 
general  wheat  yield  equation  containing  soil  water  at  seeding,  precipitation, 
and  temperature  variables  which  is  applicable  to  a  widely  representative 
soil,  optimal  decision  rules  can  be  estimated  across  the  various  locations 
and  climates  of  the   spring  wheat  region.   One  piece  of  information  needed 
to  obtain  an  estimate  of  the  decision  rule  is  average  soil  water  after  a 
year  of  fallow  (and  after  a  year  in  crop,  depending  on  the  region);  this 
information  is  substituted  into  the  wheat  yield  equation  as  soil  water 
at  seeding  to  estimate  the  right  hand  side  of  (2.16).   The  other  informa- 
tion required  for  an  estimate  of  both  sides  of  (2.16)  is  a  fairly  long 
set  of  time  series  data  on  precipitation  and  temperature  measurements 
during  the  growing  season.   A  series  as  short  as  ten  years  or  even  less 
could  be  useful,  but  30  or  40  years  would  be  aesirable. 
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Before  going  into  more  detail  on  how  the  wheat  yield  response  surface 
is  utilized  for  developing  economic  decision  rules,  the  ernpirically  esti- 
mated relationship  is  presented. 

The  Spring  Wheat  ^ield  Equatior- 
The  data  base  for  estimation  of  yield  response  was  long-term  experi- 
ments done  cooperatively  by  the  U.S.O.A.  and  state  agricultural  experi- 
ment stations  at  Mandan,  North  Dakota  and  Havre,  Montana.   The  data  for 
Montana  came  directly  from  (Army  and  Hide),  while  that  from  North  Dakota 
was  obtained  from  the  files  at  the  Mandan  Experiment  Station.   Details 
of  the  rotations  started  by  the  Office  of  Dryland  Agriculture,  U.S.D.A. 
are  given  in  (Chilcott).   The  rotations  used  for  spring  wheat  were  those 
usually  coded  MCA  (continuous  cropping)  and  MCD  (alternate  crop-fallow), 
see  for  example  (Army  and  Hide).   Soil  moisture  measurements  to  four  feet 
were  taken  on  these  treatments  by  one- foot  layers,  while  most  of  the  other 
rotations  which  might  have  been  of  interest  did  not  have  soil  moisture 
measurements  taken. 

A  great  deal  of  research  effort  was  expended  trying  to  utilize  spring 
wheat  rotation  data  from  other  locations,  but  to  no  avail.   Soil  water 
data  was  either  inadequate  or  questions  arose  as  to  the  validity  of  the 
yield  data.   At  Newell,  South  Dakota,  we  could  not  get  the  weather  data 
to  correlate  with  crop  yields  and  had  to  conclude  that  the  weather  data 
available  was  not  what  it  was  supposed  to  be. 

Even  at  Mandan,  North  Dakota,  where  we  used  the  data,  there  was  no 
record  of  emergence  dates  for  the  crops  which  reduced  the  jirccision  in 
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estimation  that  would  have  been  possible  otherwise.   We  estimated  emer- 
gence date  at  Mandan  from  the  detailed  information  available  at  the 
Havre,  Montana,  location. 

One  serious  problem  at  Havre  was  the  variation  in  soil  moisture 
measurement  dates.   Soil  water  was  supposedly  measured  at  emergence,  but 
was  actually  done  sometime  afterwards  and  this  variation  was  quite  a 
problem.   An  evaporation-precipitation  accounting  model  was  constructed 
to  adjust  measured  soil  water  to  its  equivalent  at  emergence.   Obviously, 
the  adjustment  was  subject  to  considerable  error,  but  this  was  the  best 
that  could  be  done.   Using  different  soil  water  measurement  dates  rela- 
tive to  emergence  date  was  unacceptable  because  precipitation  measured  from 
emergence  date  would  then  be  confounded  with  the  soil  water. 

At  Mandan,  soil  water  was  measured  at  seeding  time  with  considerable 
variation;  the  time  of  measurement  was  usually  after  seeding.   These 
measurements  were  adjusted  to  correspond  with  an  estimated  emergence 
date  based  on  average  time  from  seeding  to  emergence  at  Havre.   Obviously, 
the  soil  moisture  variables  which  had  to  be  used  in  the  yield  response 
analysis  were  subject  to  considerable  error,  especially  when  considered 
as  separate  variables  for  each  one-foot  layer,  as  contrasted  to  a  single 
variable  for  the  entire  four  feet. 

Evidence  which  evolved  in  fitting  separate  equations  to  the  two 

locations  suggested  that  soil  moisture  variables  were  less  dependable 

at  Mandan  than  at  Havre.   The  most  likely  explanation  for  this  difference 

in  reliability  is  the  need  at  Mandan  to  estimate  emergence  date  from 

seeding  date,  while  emergence  date  was  recorded  at  Havre. 

/ 


-63- 

The  yield  response  analysis  used  emergence  date  as  the  reference 
point  from  which  weekly  periods  were  defined  during  the  growing  season. 
Precipitation  and  temperature  variables  were  defined  for  each  of  13  weeks 
measured  from  emergence  date.   Errors  in  estimated  emergence  dates  tend 
to  confound  net  effects  of  precipitation  and  Leriiipera  1.  ure  among  the 
various  weekly  periods  of  phenological  development  of  spring  wheat.   in 
addition,  early  season  precipitation  variables  are  confounded  with  the 
soil  moisture  measurements. 

In  spite  of  these  serious  problems  with  the  data,  relatively  good 
statistical  results  were  obtained  for  the  pooled  data  at  Havre  and  Mandan, 
or  at  least,  the  statistical  results  appeared  to  be  quite  satisfactory. 
Of  course,  the  data  problems  should  not  be  ignored  in  interpreting  the 
results  because  the  consequences  of  errors  in  variables  are  not  always 
apparent  in  sample  results. 

The  algebraic  form  of  the  equation  fitted  to  the  data  is  nonlinear 
and  quite  complicated.   A  multiplicative  error  term  comparable  to  the  model 
described  in  Section  II  was  specified  for  spring  wheat.   This  form  for 
the  statistical  error  proved  superior  to  the  usual  additive  form.   One 
source  of  nonlinearity  stems  from  the  multiplicative  statistical  error 
because  the  logarithmic  transformation  required  to  get  an  additive  error 
requires  taking  the  logarithm  of  mean  yield.   Unless  mean  yield  is  mul- 
tiplicative in  the  independent  variables,  the  logarithmic  transformation 
implies  nonlinearity  in  the  parameters  of  the  model . 

A  composite  variable  is  used  which  is  a  linear  function  with  respect 
to  its  parameters  but  contains  some  interaction  terms  which  arc  nonlinear 
in  the  water  and  temperature  variables.   The  variables  are  defined  below: 
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W,  ,  W„,  W., ,  W,  are  plant  available  soil  water  variables  for  the  first 
1    2    3    4 

through  fourth  foot  of  soil; 
P, ,  P„,  .  .  .  P,^  are  weekly  precipitation  variables  ordered  from 

emergence  date; 
T-i  ,  T2,  .  •  •  Tj^3  are  weekly  average  daily  maximum  temperature  variables 

ordered  from  emergence  date; 
X;L  =  W^  +  W2  +  Pi  +  .  .  .  Pg; 

X2  =  Ty  +  Tg  +  Tg; 

X3  =  X1/X2; 

X^  =  (W^  +  ...  +  W^  +  P^  +  ...  +  P^^)^/2(^^  +  ...  +  T^2)/13 

{0  if  the  location  is  Havre 
1  if  the  location  is  Mandan 
The  X-variables  are  interactions  to  allow  a  more  complex  relationship  between 
water  and  temperature,  and  among  water  and  temperature  variables  in  different 
time  periods  within  the  growing  season.   In  an  intuitive  sense,  X-i  is  an 
early-season  water  variable;  X9  is  a  critical-development-period  temperature 
variable;  Xt  is  an  interaction  between  Xj^  and  X2;  and  X4  is  a  total  season 
interaction  between  water  and  temperature. 

The  composite  variable  which  is  specified  to  be  a  homogenous  linear 
function  of  the  above  defined  water  and  temperature  variables  is  desig- 
nated by  the  letter  u,  and  the  yield  equation  fitted  to  the  data  is 
(3.1)      y  =  aibd  -  e-^'')       ,      i  =  1,  2,  3 

where  a-,,  a2,  a-^,  b,  and  c  are  unknown  parameters.   The  three  parameters 
a,  ,  a„,  and  a-,  adjust  the  asymptotic  yield  parameter  b  according  to  whether 
the  crop  was  planted  ca   fallow  or  not  and  according  to  location  (Havre  or 
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Mandan) .   The  norm  reoresented  by  b  is  taken  as  wheat  planted  on  stubble 
at  Havre  and  tne  other  three  parameters  are  proportionate  adjustments  in 
b.   The  definitions  for  these  adjustments  in  terms  of  asymptotic  yields 
are  as  follows: 

b  corresponds  to  Havre  on  stubble 

a-,b  corresponds  to  Havre  on  fallow 

a  b  corresponds  to  Mandan  on  stubble 

a^b  corresponds  to  Mandan  on  fallow. 

The  general  shape  of  the  function  in  (3.1)  is  concave  downward  if  c 
is  less  than  one,  but  if  c  is  greater  than  one,  the  curve  has  an  inflec- 
tion point  to  the  right  of  zero.   Thus  when  c  is  greater  than  one,  the  curve 
is  the  classic  sigmoid  shape  postulated  in  production  economics. 

Each  variable  which  enters  the  composite  variable  u  linearly  also 
requires  another  parameter  to  be  estimated.   A  straightforward  count  gives 
a  total  of  40  parameters  to  represent  the  yield  equation;  5  for  (3.1) 
when  the  adjustments  to  b  are  considered  and  35  variables  listed  which 
enter  u  linearly.   Clearly  we  could  use  some  device  to  reduce  the  number 
of  parameters.   The  method  used  is  due  to  R.  A.  Fisher;  a  polynomial  con- 
straint with  time  as  the  argument  is  placed  on  the  coefficients  associated 
with  each  P,  ,  ?„,  •••  P^n  and  T,,    T^ ,  ...  Tj,^. 

Let  t  be  a  continuous  variable  associated  with  time  during  the  growing 

season  measured  from  emergence  date,  i.e. ,  t  =  0  at  emergence  date  and  t  =  13 

for  the  end  of  the  13th  week  after  emergence.   The  parameter  estimates  for 

the  linear  coefficients  on  precipitation  are  denoted  by  d  ,  d  ,  . . .  d   .   Then 

1    2       13 

these  coefficients  are  forced  to  obey  a  polynomial  as  follows: 
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(3.2)  dj.  =  Oq  +  a^t  +  a2t2    ,    t  =  1,  2,  ...  13. 
Consequently,  the  13  parameters  d^,  .  .  .  d-^2   ^^^   completely  determined 
by  the  three  parameters  a^.,  a,,  and  a^  which  reduces  the  number  of  para- 
meters in  the  model  by  ten. 

The  same  procedure  is  used  to  estimate  the  linear  coefficients  on 
temperature;  let  these  parameters  be  denoted  by  e  ,  e  ,  .  .  .  e^_.   The 
polynomial  used  to  constrain  these  parameter  estimates  is 

(3.3)  ^t  "  ^0  "^  "^1*^  "*"  ^2^^   "^  ^3*^"^     '    t  =  1,  2,  ...  13 
which  saves  nine  parameters  in  the  statistical  estimation.   A  cubic 
instead  of  a  quadratic  was  required  on  temperature  variables  to  provide 
adequate  flexibility. 

In  total,  the  polynomial  constraint  technique  saved  19  parameters, 
giving  a  total  net  number  of  parameters  in  the  model  equal  to  21.   After 
estimation  of  the  polynomial  parameters  in  (3.2)  and  (3.3),  the  individual 
weekly  coefficients  for  precipitation  and  temperature  can  be  calculated 
as  well  as  statistical  standard  errors  for  them  because  the  transforma- 
tions involved  are  linear. 

Numerical  results  for  the  nonlinear  regression  estimation  are  given 
in  Table  14.   The  asymptotic  yield  parameter,  denoted  by  b  in  (3.1),  tended 
to  increase  without  limit  in  the  statistical  fitting  procedure.   There- 
fore this  parameter  was  arbitrarily  set  equal  to  40  waich  gave  about  the 
same  fit  as  very  large  numbers  obtained  from  the  fitting  procedure. 
Apparently  the  data  do  not  contain  high  enough  yields  to  estimate  this 
parameter. 
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Table  14 


Statistical  Results  for  Spring  Wheat  Yield  Response 


Description  of  Parameter 
or  Variable 


Regression 
Coefficient 


Standard 

Error 


Asymptotic  yield 

Adjustments  in  asymptotic  yield 

Fallow  at  Havre 

Stubble  at  Mandan 

Fallow  at  Mandan 


40 


t-ratio 


1.22 

0.0721 

3.04* 

1.09 

0.0892 

1.06* 

1.30 

0.0928 

3.21* 

Exponent  on  composite  variable  u 


0.973 


0.114 


0.24* 


Soil  Water  Variables 
1st  foot 
2nd  foot 
3rd  foot 
4th  foot 

Precipitation  Variables 

1st  week 

2nd  week 

3rd  week 

4th  week 

5th  week 

6th  week 

7th  week 

8th  week 

9th  week 
10th  week 
11th  week 
12th  week 
13th  week 

Temperature  Variables 

1st  week 

2nd  week 

3rd  week 

4th  week 

5th  week 

6th  week 

7th  week 

8th  week 

9th  week 
10th  week 
11th  week 
12th  week 
13th  week 


145 

.0473 

3.07 

185 

.0422 

4.38 

129 

.0385 

3.35 

114 

.0441 

2.59 

168 

.0453 

3.71 

164 

.0433 

3.79 

158 

.0414 

3.82 

151 

.0394 

3.83 

14  3 

.0374 

3.82 

133 

.0354 

3.76 

123 

.0334 

3.68 

112 

.0315 

3.56 

100 

.0300 

3.33 

086 

.0290 

2.98 

072 

.0288 

2.49 

056 

.0297 

1.89 

040 

.0318 

1.25 

00161 

.00111 

-1.45 

0004  2 

.00080 

-0.53 

0002  7 

.0006  7 

0.40 

00054 

.00059 

0.92 

0004  7 

.0005  3 

0.89 

00014 

.00052 

0.27 

00035 

.00059 

-0.59 

00094 

. 00069 

-1.36 

0015  3 

.00076 

-2.01 

00205 

.00074 

-2.77 

00241 

.00063 

-3.83 

00253 

.00063 

-4.02 

00233 

.00113 

-2.06 
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Table  14  (cont'd) 


Description  of  Parameter 
or  Variable 


Regression 
Coefficient 


Standard 

Error 


t-ratio 


Interaction  Variables 
Early  season  water,  X 
Critical  season  temperature,  X2 
X]^/X2 
Total  water-temperature,  X^ 

Constant  for  Mandan  location,  Z 


.357 

.0773 

-4.62 

.0194 

.00468 

4.15 

29.4 

6.24 

4.71 

.00658 

.00198 

-3.32 

.0903 


.0276 


3.27 


*The  null  hypothesis  associated  with  the  t-statistic  is  that  the  parameter 
is  equal  to  one  (instead  of  zero  in  the  other  cases). 
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Ap  a  point  of  reference  in  the  phenology  cf  the  crop,  heading  is 
about  the  middle  of  week  number  eight  and  the  soft-dough  stage  is  about 
the  middle  of  week  number  twelve.   However,  the  coefficients  on  the 
weekly  precipitation  and  temperature  variables  cannot  be  interpreted  in 
isolation  from  the  interaction  terms  where  these  same  variables  appear. 

Net  effects  of  individual  variables  (partii'l  derivatives)  are  given 
in  Table  lb    for  all  variables  at  their  mean  levels  for  the  entire  sample 
of  data  at  both  locations.   The  quite  large  detrimental  effects  for 
temperature  during  weeks  7,  8,  and  9  stem  from  the  choice  of  interaction 
terms,  namely,  the  ratio  of  early  season  water  to  critical  season  tempera- 
ture.  Choice  of  the  critical  season  for  temperature  was  based  on  a 
simple  comparison  of  periods  around  heading  time,  and  there  is  no  reason 
to  expect  that  the  net  effects  are  actually  as  choppy  as  estimated  in  Table 
15.   Nevertheless,  the  net  effects  for  temperature  would  appear  quite 
plausible. 

One  might  have  expected  net  effects  of  precipitation  to  be  rela- 
tively larger  in  weeks  6,  7,  and  8  preceding  heading,  but  apparently 
early  season  precipitation  is  the  most  important  when  recognition  is 
made  of  the  empirical  distribution  of  precipitation  for  the  region  under 
study.   The  nearly  zero  and  negative  effects  for  precipitation  in  weeks 
11  through  13  suggest  that  the  crop  cannot  effectively  use  water  that  late 
in  the  season.   However,  these  effects  are  increased  by  the  total  season 
water-termperature  interaction  when  total  water  is  relatively  large; 
possibly  the  extra  water  delays  maturity  and  thus  makes  water  later  in 
the  season  more  productive. 
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Table  15 


Individual  Variable  Net  Effects  on  Yield  With 
All  Variables  at  Mean  Levels 


Description  of  Variable 


Net  Effect  on  Yield* 
(bushels  per  acre) 


Soil  Water 
1st  foot 
2nd  foot 
3rd  foot 
4th  foot 


1.87 
2.71 
1.54 
1.23 


Precipitation  Variables 

1st  week 

2nd  week 

3rd  week 

4th  week 

5th  week 

6th  week 

7th  week 

8th  week 

9th  week 
10th  week 
11th  week 
12th  week 
13th  week 


2.35 
2.27 
2.14 


00 
83 
62 
42 
19 


0.56 

0.27 

-0.02 

-0.35 

-0.69 


Temperature  Variables 

1st  week 

2nd  week 

3rd  week 

4th  week 

5th  week 

6th  week 

7th  week 

8th  week 

9th  week 
10th  week 
11th  week 
12th  week 
13th  week 


-0.50 
-0.48 
-0.46 
-0.46 
-0.46 
-0.47 
-1.08 
-1.09 
-1.10 
-0.51 
-0.52 
-0.52 
-0.52 


*Net  effect  on  yield  is  technically  the  partial  derivative  of  the  yield 
equation  with  respect  to  the  designated  variable  with  all  variables  at 
their  sample  means. 
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Net  effects  for  soil  water  by  one-foot  layers  seems  to  make  sense 
on  a  priori  reasoning.   The  second  foot  has  a  relatively  high  density 
of  roots  and  is  protected  from  evaporation  losses  compared  to  the  top 
foot.   The  3rd  and  4th  feet  are  progressively  less  effective  sources 
of  water  because  of  fewer  roots  and  greater  distance  from  the  main 
plant  body. 

Supplementary  statistical  results  are  given  In  Table  16.   The 
statistical  error  correlations  among  treatments  are  recognized  in  the 
statistical  model  and  the  reported  correlations  are  maximum  likelihood 
estimates.   Note  the  high  correlations  between  stubble  and  fallow 
plots  at  the  same  locations  as  would  be  expected.   The  model  seemed  to 
fit  the  fallow  treatment  at  Havre  the  worst,  as  evidenced  by  the  larger 
standard  error  of  estimate  and  smaller  Durbin-Watson  statistic;  the 
latter  measures  serial  correlation  in  the  time  series  of  residuals  from 
the  regression  equation  for  the  specific  treatment. 

Estimation  of  Conditional  Expected 
Yields  in  Relation  to  Soil  Water 

The  empirical  estimate  of  the  Havre-Mandan  yield  equation  given 
in  Table  14  can  be  used  together  with  a  time  series  on  climatic  data  at 
a  given  site  to  estimate  expected  yield  of  spring  wheat  as  a  function 
of  plant  available  soil  water  at  seeding  time.   This  functional  relation- 
ship will  be  referred  to  as  simply  "the  mean  yield  function  in  soil 
moisture. " 
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Table  16 
Supplementary  Statistical  Results  for  Spring  Wheat 


Standard  Error  of  Estimate* 
Durbin-Watson  Statistic 
Statistical  Error  Correlations 

Havre  Stubble 

Havre  Fallow 

Mandan  Stubble 

Mandan  Fallow 


Treatments 

Havre 

Havre 

Mandan 

Mandan 

Stubble 

Fallow 

Stubble 

Fallow 

.241 

.304 

.220 

.237 

2.41 

1.39 

2.35 

2.22 

1.0 

.571 

-  .041 

-  .114 

1.0 

.275 

-  .055 

1.0 

.427 
1.0 

Degrees  of  Freedom  =  98 

Coefficient  of  Multiple  Determination  (R^)  =  .926 


*Unit  of  measure  is  the  natural  logarithm  of  yield.   Comparable  measures 
in  the  original  units  of  bushels  per  acre  for  Havre  and  Mandan  are  3.26 
and  3.80,  respectively,  at  mean  yields  for  each  location.   (The  average 
standard  errors  across  stubble  and  fallow  were  used  to  get  these  esti- 
mates. ) 
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This  functional  relationship  could  be  estimated  with  four  variables, 
one  each  for  the  four  layeis  of  soil  water  by  one-foot  depths,  but  the 
economic  analysis  was  not  attempted  with  that  much  precision  on  soil 
water.   The  primary  reason  for  not  doing  so  was  its  impracticality  in 
application  to  actual  farm  decisions;  farmers  would  probably  use  a  soil 
probe  to  estimate  total  water  and  tne  precision  of  such  a  procedure  would 
not  warrant  differentiation  of  water  among  the  one-foot  soil  layers. 
In  addition,  the  economic  analysis  is  much  simpler  and  more  easily  pre- 
sented if  there  is  only  one  soil  moisture  variable. 

At  Havre  and  Mandan,  data  were  available  by  one-foot  depths  for  an 
extended  number  of  years.   Mean  plant  available  water  in  each  foot  of 
soil  was  used  as  a  relative  weight  to  reduce  the  individual  soil  water 
regression  coefficients  in  Table  14  to  a  single  coefficient  for  soil 
water  to  four  feet.   The  final  coefficient  is  a  weighted  average  of  the 
one-foot  coefficients,  with  mean  soil  water  by  one-foot  depths  serving  as 
weights  at  a  given  location.   At  Dickenson,  North  Dakota  and  Sidney, 
Montana,  the  weighted  coefficient  at  Mandan  was  used  in  the  analysis.   It 
was  obvious  from  the  Havre  and  Mandan  work  that  the  weighted  average 
coefficient  was  not  very  sensitive  to  different  locations. 

Suppose  that  n  years  of  climatic  data  are  available  at  a  given  site, 
i.e.,  daily  maximum  temperature  and  precipitation.   A  specific  emergence 
date  is  specified  for  the  location  from  which  13  weekly  temperature  and 
precipitation  variables  are  constructed  for  each  of  the  n  years  of  data. 
These  weekly  climatic  variables  are  used  to  construct  the  interaction 
variables  in  the  spring  wheat  yield  equation  (see  tlie  definition  of 
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variables  just  before  (3.1)). 

The  next  step  in  calculating  the  mean  yield  function  in  soil  moisture 
is  to  estimate  the  most  appropriate  asymptotic  yield,  a.b,  to  use  in  (3.1). 
The  only  choice  at  Mandan  or  Havre  is  between  planting  on  stubble  or  fallow. 
The  assumption  was  made  that  the  fallow  asymptotic  yield  parameter  is 
most  appropriate  for  planting  on  either  fallow  or  stubble  when  nitrogen 
fertilization  is  practiced;  adjustments  are  made  in  production  costs  for 
the  difference  in  nitrogen  applied  to  stubble  versus  fallow.   The  statistical 
work  for  the  Highwood  Bench  data  in  Section  II  suggested  that  there  exists 
a  "pure  fallow  effect"  in  addition  to  the  influence  of  soil  water  and  plant 
available  nitrogen,  but  we  cannot  separate  such  an  effect  out  from  the 
nitrogen  effect  in  these  data  from  Havre  and  Mandan  because  nitrogen 
fertilization  was  not  a  factor  in  the  experiments.   At  Dickenson,  North 
Dakota  and  Sidney,  Montana,  the  fallow  asymptotic  yield  parameter  at 
Mandan  was  assumed  to  be  applicable.   It  is  noted  that  the  difference  in 
asymptotic  yield  estimated  for  planting  on  fallow  at  Mandan  is  only  seven 
percent  higher  than  at  Havre,  certainly  not  a  very  significant  consideration. 

The  primary  calculations  to  estimate  the  mean  yield  function  in  soil 
moisture  are  as  follows.   For  a  given  level  of  plant  available  soil  moisture, 
the  yield  is  calculated  from  the  estimated  parameters  of  (3.1)  as  given  in 
Table  14  (with  selections  and  modifications  described  above)  for  each  of 
the  n  years  of  climatic  data  available.   Let  us  denote  a  given  level  of 
plant  available  soil  water  in  the  top  four  feet  of  soil  by  W.   Then  for  a 
given  value  of  W,  wc  have  calaculated  n  predicted  spring  wheat  yields,  one 
for  each  year  of  the  time  series.   The  simple  average  of  those  n  yields 
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is  the  estimated  conditional  mean  yield,  given  soil  water  at  planting  time* 
equals  W   The  standard  deviation  of  conditional  mean  yields,  given  W,  can 
be  calculated  from  the  same  n  predicted  yields  by  the  usual  formula  from 
statistical  theory. 

These  conditional  mean  yields  and  standard  deviations  are  given  in 
tabular  form  in  the  results  reported  below  for  each  location.   For  example, 
the  Havre  results  are  given  in  Table  17.   Conceptually,  conditional  mean 
yields  and  the  standard  deviation  of  yields  are  each  a  continuous  function 
of  plant  available  water  at  seeding  time.   If  desired,  a  continuous  function 
could  be  fitted  to  the  tabular  data  by  regression  to  get  an  estimate  of  the 
corresponding  continuous  functions,  but  the  tabular  approximations  are 
adequate  for  our  purposes. 

Cropping  Strategy  Decision  Models 
For  Spring  Wheat  Areas 

Two  types  of  decision  models  were  used  to  determine  optimal  cropping 

strategies  for  spring  wheat  in  the  Northern  Great  Plains.   The  more  complex 

model  is  a  mathematical  model  known  as  a  Markov  chain.   A  simpler,  more 

intuitive  model  used  was  the  "basic  trade-off  equation"  (BTE)  described 

earlier  in  the  text. 

Markov  Chain  Decision  Model 

The  basic  decision  model  to  analyze  the  choice  between  whether  to 
plant  a  spring  grain  crop  or  to  fallow  the  land  is  mathematically  de- 
scribed as  a  Markov  chain.   A  complete  explanation  of  the  model  in  the 
context  used  here  is  given  in  (Burt  and  Allison,  Burt  and  Stauber). 


*Plant  available  soil  water  at  emergency  is  assumed  to  be  approximately 
equal  to  the  amount  measured  at  seeding  time. 
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Briefly,  the  problem  is  to  make  a  decision  concerning  the  use  of 
land  (crop  or  fallow)  that  will  maximize  present  value  of  net  returns  per 
acre  over  a  given  planning  horizon.   In  the  model  the  decision  regarding 
the  utilization  of  a  tract  of  land  is  made  shortly  before  seeding  time  for 
spring  wheat.   At  this  decision  point,  a  quantitative  evaluation  is  made 
of  the  information  concerning  plant  available  moisture.   This  evaluation 
is  used  as  the  basis  for  making  the  decision  on  land  use.   That  is,  infor- 
mation regarding  plant  available  moisture  at  seeding  time  is  used  as  the 
basis  for  deciding  whether  to  plant  spring  wheat  or  to  leave  the  land 
fallow. 

In  our  present  situation  we  are  concerned  only  with  a  model  to  deter- 
mine whether  the  land  should  be  fallowed  or  planted  to  spring  wheat.   Assum- 
ing that  consecutive  years  of  fallow  are  not  permitted  for  a  given  tract 
of  ground,  a  cropping  strategy  would  be  comprised  of  the  following  simple 
rule: 

If  the  land  was  in  fallow  last  year,  plant  wheat  this  year. 
If  the  land  was  in  crop  last  year,  determine  the  amount  of 
available  soil  water  at  planting  time.   If  available  soil 
water  is  M^  inches  or  greater  plant  wheat,  but  if  available 
soil  water  is  less  than  Mg  inches,  fallow  the  land. 

The  critical  available  soil  moisture  value  M  is  simply  a  given 
number  of  inches  of  plant  available  water  and  it  completely  determ.ines 
the  decision  rule.   Conceptually,  there  is  some  optimal  choice  of  M 
for  a  given  location,  price  of  wheat,  costs  of  production,  and  technology 
of  production. 

In  order  to  find  the  optimal  value  of  i^. ,  in  a  given  situation,  certain 
types  of  information  must  be  available  for  specii icd  levels  of  available  soil 
moisture.   Therefore,  tlie  possible  values  of  plant  available  moisture  are 
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defined  as  discrete  intervals  which  are  designated  by  the  midpoints  of  the 
intervals.   For  instance,  plane  available  soil  moisture  between  the  values  of 
2.0  inches  to  3.0  inches  is  designated  by  the  midpoint  of  the  interval, 
2.5  inches.   In  the  Markov  chain  model,  a  variable  such  as  soil  moisture 
which  describes  the  condition  of  a  decision  process  is  called  a  state  variable. 

Two  general  types  of  empirical  information  are  needed  to  estimate 
reliable  cropping  strategies  for  a  given  location.   One  necessary  empirical 
relationship  is  an  estimate  of  average  crop  yield  in  relation  to  available 
soil  moisture  at  seeding  time.   That  is,  for  a  given  crop,  we  need  to  know 
average  yield  over  many  years  of  climatic  variation  for  a  given  level  of 
available  soil  moisture  at  seeding  time.   The  second  type  of  empirical 
relationship  needed  is  the  probability  distribution  for  available  soil 
moisture  changes  from  one  season  to  the  next  under  given  land  use  alter- 
natives.  If  both  types  of  empirical  information  are  available,  the  more 
complex  Markov  chain  decision  model  can  be  applied.   If  the  conditional 
expected  yield  relationship  is  available  the  basic  trade-off  equation  can 
be  applied  even  though  the  required  probability  distribution  on  soil 
moisture  is  not  available. 

The  necessary  data  were  available  to  apply  the  Markov  chain  model  at 
Havre,  Montana,  and  Mandan,  North  Dakota.   Probability  distributions  were 
not  available  at  Sidney,  Montana  and  Dickinson,  North  Dakota,  so  the  BTE 
was  used  at  these  locations.   The  BTE  was  used  at  Havre  and  Mandan  also 
to  provide  a  comparison  of  results  between  the  two  decision  models. 

A  summary  of  the  empirical  data  used  in  the  models  and  the  economic 
results  of  the  models  are  summarized  below  by  locations. 
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Havre,  Montana 

The  expected  yield  relationship  for  Havre  was  estimated  by  using  the 
Havre-Mandan  yield  equation  reported  in  Table  14  together  with  a  series 
of  weather  data  from  Havre.   The  conditional  yield  estimates  and  their 
standard  deviations  are  presented  in  Table  17. 

Probability  distributions  for  available  soil  moisture  transitions 
from  one  season  to  the  next  under  given  land  use  alternatives  were 
obtained  from  recorded  soil  moistures  at  Havre.   A  formal  statement  of 
the  model  used  to  relate  soil  moisture  at  period  t  to  soil  moisture  at 
period  t-1  is  given  by 
(3.4)      Mj.  =  f(M^_j^,  LUj__^,  LU^.) 

Literally  this  relationship  states  that  available  soil  moisture  at 
the  end  of  the  present  time  period  (Mj.)  is  a  function  of  soil  moisture  at 
the  end  of  the  previous  time  period  (Mj-.j),  land  use  during  the  previous 
time  period  iUJ^_i)    and  land  use  during  the  present  time  period  (LU  ). 
Where  crop  (C)  or  fallow  (F)  are  the  land  use  alternatives,  four  situations 
exist:  they  can  be  enumerated  as: 

1.  M,  =  f(M,_i,  F^_^.  F^) 

2.  M^_  =  f(M^_i,  F^_^,  C^) 

3-    \   =  f(Vl'  Vl'  ^t^ 

'■       \   =  '^\-l'   Vl'  ^) 

Linear  regression  techniques  were  used  to  estimate  the  relationship 
between  M^.  and  M|._^  for  situations  2,  3,  and  4.   No  data  were  available  to 
estinate  the  relationship  depicted  in  situation  1.   It  was  assumed  that 
the  estimates  obtained  for  situation  ^  would  suffice  for  situation  1.   This 
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Table  .17 


Conditional  Expected  Yields  of  Spring  Wheat  for  Specified  Levels 
Of  Available  Soil  Moisture  at  Havre,  Montana 


Available  Conditional  Standard 

Soil  Moisture  Expected  Yield        Deviation 


(inches)  (bu)  (bu) 

.5  3.98  5.99 

1.0  4.83  6.44 

1.5  5.89  6.88 

2.0  7.05  7.22 

2.5  8.31  7.46 

3.0  9.68  7.56 

3.5  11.07  7.63 

4.0  12.45  7.67 

4.5  13.82  7.68 

5.0  15.17  7.66 

5.5  16.49  7.63 

6.0  17.78  7.57 

6.5  19.04  7.51 

7.0  20.27  7.42 

7.5  21.47  7.32 

8.0  22.62  7.20 
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assumption  depends  on  no  significant  difference  in  the  distribution  of 
M    throughout  the  soil  profile  following  fallow  and  crop.   The  validity 
of  this  assumption  is  questionable  but  lack  of  data  at  Havre  or  similar 
sites  left  no  viable  alternative. 

The  regression  model  to  be  estimated  for  a  given  situation  is  expressed 
as 
(3.5)     Mj.  =  a  +  g  M^_^  +  e^ 

Regression  estimates  for  situations  where  data  were  available  are 
summarized  in  Table  18. 

Table  18 

Estimates  of  Parameters  for  Regression  of  Available 

Soil  Water  in  Time  t  on  Available  Soil  Water 

In  Time  t-1  at  Havre,  Montana 


Land  Use  Estimates 

Situation  ^  g 


FC  4.025  .3650  1.320 

CC  5.608  .1271  1.150 

CF  7.335  .2311  1.123 

The  immediately  following  discussion  treats  the  estimated  unknown 
parameters  in  the  regression  model  as  though  they  were  the  actual  values 
because  the  decision  model  takes  these  sample  estimates  as  the  "true"  values. 
Under  the  assumptions  of  the  linear  regression  model,  Mj.  is  normally 
distributed  witli  mean,  E(Mj.)  =  a  +  S^t-i  and  variance,  Var  (Mj.)  = 
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E(e{-  )  =  a^.  This  information  is  used  to  derive  the  conditional  prob- 
abilities needed  in  the  Markov  cliain  decision  model.  7'he  procedure  is 
sketched  briefly  h'Te. 

A  given  value  ol'  Mj-_]  implies  the  mean  value  of  M^  ,  i.e.,  l^fM^) 
a   +  3Mj._T  .   Since  M^.  is  distributed  normally  with  standard  deviation,  o, 

the  probability  of  occurrence  of  any  specified  value  for  Kj-  can  be  found 
using  the  Z  transformation, 

Mt  -  E(M^) 
—O =^ 

and  the  standard  normal  tables. 

Cost  and  Returns  Data 

The  cost  data  for  the  land  use  alternatives,  summer  fallow  or  planting 
spring  wheat^ were  taken  from  cost  and  returns  studies  for  McCone  and  Sheridan 
counties  (Schaefer,  Stauber,  and  Luft,  1976a,  1976b). 

The  variable  costs  were  adjusted  to  reflect  lower  levels  of  input  use 
at  unfavorable  soil  moisture  levels.   The  details  of  the  adjustments  are 
not  presented  but  in  general  the  input  reductions  were  assumed  to  occur  in 
seed,  nitrogen  and  phosphorus  fertilizer,  insecticide,  miscellaneous  expense 
and  interest  on  operating  capital.   Crop  insurance  was  omitted  as  a  cost 
item  because  the  P'CIC  does  not  issue  loans  on  recrop  acreages.   Charges  for 
operator  labor  were  also  omitted  from  the  variable  costs.   The  costs  used  at 
Havre  are  presented  in  Table  19. 

The  greater  variable  cost  of  spring  wheat  on  stubble  compared  to  fallow 
is  mainly  due  to  extra  nitrogen  fertilizer  costs  to  compensate  for  lack 
of  mineralized  nitrogen  under  recrop  conditions. 
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Table  19 


Yield  and  Variable  Costs  by  Specified  Soil  Moisture  Levels 
For  Recrop  Spring  Wheat,  Spring  Wheat  on  Fallow 
And  Fallow  at  Havre,  Montana 


Soil 
Moisture 


Yield 


Variable  Costs  of 
Spring  Wheat 


Fallow 


Recrop 


Variable  Costs 
of  Fallow 


(inches) 

(bushels) 

(dollars) 

(dollars) 

(d( 

ollars) 

.5 

3.98 

5.00 

5.00 

4.00 

1.5 

5.89 

10.00 

10.00 

4.00 

2.5 

8.31 

15.00 

16.70 

4.00 

3.5 

11.07 

20.00 

22.25 

4.00 

4.5 

13.82 

20.00 

25.25 

4.00 

5.5 

16.49 

20.00 

25.25 

4.00 

6.5 

19.04 

20.00 

25.25 

4.00 

7.5 

21.47 

20.00 

25.25 

4.00 
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Gross  returns  and  expected  Litunediate  returns  per  acre  were  calculated 
for  three  situations  as  enumerated  below  and  presented  in  Table  20: 

1.  A  wheat  price  of  $3.00  per  bushel  and  100  percent  of  estimated 
yields 

2.  A  wheat  price  of  $3.00  per  bushel  and  120  percent  of  estimated 
yields 

3.  A  wheat  price  of  $3.50  per  bushel  and  120  percent  of  estimated 
yields 

Economic  Results 

The  results  of  the  Markov  Chain  Model  are  presented  in  Table  21 
and  graphically  in  Figure  1.   The  numerical  results  were  derived  by  solving 
the  system  of  simultaneous  equations  (Burt  and  Allison,  p.  128)  implied 
by  selected  values  of  M  ,  the  critical  value  for  available  soil  moisture. 
Recall,  that  if  available  soil  moisture  is  M  inches  or  greater,  wheat 
is  planted.   If  available  soil  moisture  is  less  than  M  ,  the  land  is 
fallowed.   For  example,  M^  =  0  implies  continuous  cropping  while  M  =  «> 
implies  an  alternate  crop-fallow  rotation,  if  consecutive  fallow  is 
excluded  on  a  given  tract  of  land. 

The  net  returns  figures  presented  in  Table  21  and  Figure  1  are, 
more  precisely,  the  amortized  expected  present  value  of  net  returns 
associated  with  an  infinite  planning  horizon  using  an  interest  rate  of 
10  percent.   The  returns  represent  the  estimated  annual  net  returns 
above  variable  costs  of  production  at  the  specified  yield  and  price 
assumptions. 

The  optimal  value  of  the  decision  variable  Mg  is  seen  to  be  4  for 


d 

n) 

4J 

c 

^] 

u 

O      " 

IM      0) 

vj 

CO    > 

c   « 

H   BC 

3 

•u  u 

<U     Cfl 

(iJ 

-d 

0)     0) 

■u    a 

nj    D. 

o 

•H    o 

(N 

XI    u 

<u   o 

d) 

B     0) 

.H 

B  a: 

J3 

i-i 

CO 

XI 

H 

T3      C 

ai    CD 

4-1 

o    S 

OJ     0 

D.  M 

X    .-( 

W     cfl 

tn 

•a 

c  c 

<fl    o 

tn   -u 

to    CO 

o   ai 

o  3 

M 

c 

■r-l 

i-i 

a. 

en 

CO 

c 
o 

XI 

0)  to 
JJ  c 
a   u 

&  m 

tu 
•u 

CO 
•H 
T3 

0) 


, 

tn 

en 

I-I 

<«• 

tl 

1 

t— 1 

6^ 

.-1 

O 

o 

CN 

■o 

i-H 

O 

O 

^^ 

• 

tn 

(^ 

I-I 

<y> 

cfl 

1 

r-l 

B-S 

.-1 

O 

O 

CM 

X) 

.H 

^w' 

O 

O 

-•— S 

. 

to 

<n 

I-I 

<«■ 

cfl 

1 

M 

^s 

iH 

o 

O 

o 

XI 

O 

LTl 

--^ 

• 

tn 

m 

M 

</> 

to 

1 

1— I 

fr« 

rH 

O 

O 

CN 

XI 

iH 

^^ 

O 

O 

."^ 

• 

to 

CO 

U 

<l> 

to 

1 

1— 1 

s^ 

r-H 

o 

O 

CM 

X! 

^ 

o 

o 

y^ 

• 

to 

n 

V-i 

o 

to 

1 

i-H 

s^ 

r-l 

o 

O 

o 

XI 

X3 

rH 

iH 

OJ 

0) 

x: 

•H 

tn 

>< 

3 

tu 

iH 

4) 

,-^ 

X> 

)-i 

tn 

to 

iH 

3 

41 

i-H 

•H 

4-> 

J= 

•H 

O 

to 

t) 

to 

C/3 

•T-l 

C 

> 

o 

•H 

< 

s 

^^ 

CM<fO-3-CTi.HCNtN 

f^r^cNeNr^Or~.tj\ 


-a-  00  -3-  CM  -^r  ~3-  v* 

rH  ,— I  CM  CO  -<f  in  vD 


en  O  CN  O  O 

en  CM  CM  ^  Lo 


CM  O 


t3\  ,H  tn 


r^  -a-  ^  m  CM 
I-H  CM  en  <■  u-1 


>3-r-.envOi-(CMi^vO 

CTi\OCMCT\CMCN00i— I 


vOr--tx)OvO-a-i— itT\ 
.-I  r-H  CM  m  m 


o 

in 

tn 

m 

u 

CM 

~d- 

o 

<J> 

<3- 

vO 

r^ 

r^ 

<y> 

to 

r^ 

r-~ 

CJ\ 

^ 

o 

CM 

o\ 

rH 

S-5 

rH 

iH 

-vT 

ON 

VC 

00 

(Ti 

c^ 

cj 

O 

0 

fH 

iH 

.H 

CM 

cn 

~a- 

in 

r~~ 

CM 

XI 

iH 

O 

O 

to 

en 

u 

m 

O 

CM 

in 

in 

vO 

<r 

e3\ 

1 

cfl 
.H 

en 

CM 

CTi 

c» 

r^ 

r^ 

in 

CM 

^s 

iH 

CT\ 

•-{ 

~a- 

CTv 

ON 

o\ 

00 

r^ 

o 

O 

.H 

■H 

rH 

CSJ 

cn 

~* 

in 

CN 

X3 

iH 

v^ 

O 

O 

-^ 

• 

tn 

en 

I-I 

-* 

r^ 

en 

rH 

vO 

r-~ 

CM 

rH 

1 

to 

tH 

CJv 

\D 

CJN 

CM 

~a- 

~3- 

rH 

-a- 

s-s 

rH 

vD 

r^ 

tJ^ 

cn 

iH 

a\ 

r^ 

-.T 

o 

o 

rH 

CM 

c^ 

cn 

<r 

o 

T3 

cM-3-ocj>-a-vDr-~r^ 

r-~P~-CTi-a-OCM<jNrH 


vO~3-s3-vDOOO\0\0 

rHCMcn~a-inv£)r^cy> 


cnocMininvD^ON 
encMo>oor~-cnincM 


~3-rHC7NO\CT\cyiCOr-~ 

rHCMCNcn-a-invor^ 


>3-r-»tnrH\or^cMrH 
erivOcjNCM-*~d-rH»3- 


rHr^<j-cnrHtyir^--a- 
iHrHCMcn-cr-j-invD 


ooo\rHr^cMCTN<-r^ 
cy.oocnooo<-o<- 


cnmoOrHcnvDONrH 

rH    rH    rH    rH    CM 


inininininLninin 
iH  CM  m  •a  in  vo  r~- 


-85- 


Table  21 


Average  Annual  Net  Returns  in  Relation  to  the 
Decision  Rule  Variable,  Mg  for  Three  Specified  Yield-Price 
Situations,  Havre,  Montana 


Decision 

Rule 
'  "s 

Average 

Annual  Net 

Ret 

urns** 

Variable 

100%-$3.00 

120%-$3.00 

120%-$3.50 

(inches  available 

soil  water) 

(dollars) 

0 

10.31 

15.84 

21.38 

1 

10.  74 

16.51 

22.29 

2 

11.52 

17.51 

23.49 

3 

12.49 

18.43 

24.40* 

4 

12.97* 

18.66* 

24.35 

5 

12.84 

18.22 

23.60 

6 

12.42 

17.57 

22.73 

7 

12.19 

17.26 

22.32 

8 

12.13 

17.17 

22.21 

*Optimum  net  returns  per  acre. 

**A  10  percent  interest  rate  was  used  for  amortization. 
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Figure  1   Average  Annual  Net  Returns  in  Relation  to  the 
Decision  Variable,  M. ,  at  Havre,  Montana 
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both  of  the  situations  where  price  of  wheat  is  equal  to  $3.00.   The  opti- 
mal value  of  the  decision  variable  is  Mg  =  3  for  the  120%  yield-$3.50  price 
situation.   In  all  cases,  alternate  crop-fallow  (M^_  =  «>)  is  economically 
superior  to  continuous  cropping  (Mg  =  0).   The  superiority  of  alternate 
crop-fallow  over  continuous  cropping  is  inversely  related  to  the  magnitude 
of  net  returns  as  affected  by  yield  and  price  assumptions. 

In  terms  of  saline  seep  abatement,  it  can  be  seen  that  the  grain  grower 
could  move  to  a  strategy  which  would  result  in  less  frequent  fallow 
without  a  great  loss  in  economic  returns.   In  the  100%  yield-$3.00  price 
situation,  fallowing  if  available  soil  water  is  3  inches  or  less  rather 
than  at  4  inches  would  reduce  annual  average  net  returns  by  $.48  per  acre. 
Lowering  the  decision  rule  value  of  M   from  3  inches  to  2  inches  however  would 
result  in  an  additional  loss  of  $.97  per  acre.   It  should  be  noted  that  this 
loss  would  occur  on  all  crop  acres  in  the  farming  unit  and  is  not  a  trivial 
loss.   A  loss  of  $.97  per  acre  is  approximately  7.5  percent  of  the  net  return 
under  an  optimal  selection  of  M  .   (The  percentage  loss  is  substantially 
more  if  fixed  costs  are  taken  into  account.  )   The  farm  operator  must  weigh 
the  direct  economic  losses  resulting  from  deviation  from  the  optimal 
decision  rule  value  for  M  against  the  potential  benefits  realized  from 
saline  seep  control. 

If  the  decision  maker  is  willing  to  sacrifice  economic  returns,  in 

order  to  control  saline  seep,  adjustment  of  the  value  of  the  decision 

variable  M   to  3  inches  has  obvious  advantages.   This  adjustment  will  cost 
s 

him  $.48  per  acre  in  net  returns  but  it  will  result  in  keeping  the  soil 

in  a  considerably  drier  condition.   It  can  be  shown  that  the  expected  amount 
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of  available  water  in  the  soil  profile  prior  to  a  year  of  fallow  would  be 
reduced  by  approximately  .6  of  an  inch  by  adjusting  the  value  of  M  down- 
ward from  4  inches  to  3  inches.   In  order  to  illustrate  the  advantage  of 
the  decision  rule  (Mg  =  3)  over  alternate  crop-fallow,  the  probability 
of  fallowing  when  available  soil  water  exceeds  3  inches  in  the  crop- 
fallow  rotation  must  be  known.   The  probability  distribution  of  soil 
water  after  a  year  of  crop  when  following  a  rotation  of  alternate  crop- 
fallow  is  given  in  Table  22. 

The  level  of  M   used  in  the  decision  rule  is  the  upper  limit  on  soil 
s 

water  at  which  fallow  is  allowed.   That  is,  if  the  rule  Mg  =  3  is  followed 
the  land  will  never  be  fallowed  when  available  soil  moisture  exceeds  3 
inches  at  seeding  time  for  spring  wheat.   Examination  of  Table  22  reveals 
the  saline  seep  abatement  advantages  of  the  rule  M  =  3.   Under  an  alternate 
crop-fallow  system  the  soil  would  be  fallowed  when  available  soil  moisture 
was  greater  than  3  inches  64  percent  of  the  time.   The  reduction  in  deep 
percolation  losses  of  soil  water  under  the  decision  rule  M  =  3j  as  con- 
trasted to  alternate  crop-fallow,  is  apparent. 

It  can  be  shown  that  the  decision  rule  Mg  =  3  will  result  in  the  farm 
being  in  crop  approximately  two- thirds  of  the  time  and  in  fallow  one-third 
of  the  tim.e.   Under  this  rule,  tne  land  is  recropped  approximately  one-half 
of  the  time.   That  is,  approximately  one-half  of  the  years  after  a  crop  is 
grown  the  land  will  be  recropped  and  i-pproximately  one-li.ilf  of  the  years  after 
a  crop  is  grown  the  land  will  he  falJowed.   Ui'der  the  assumptions  of  the 
model ,  a  crop  will  always  be  planted  after  a  year  of  fallow. 

It  is  important  to  recognize  that  not  only  is  the  frequency  of  crop 
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Table  22 


Probabi-lity  Distribution  of  Available  Soil  Water 
After  a  Year  of  Crop  for  Alternate 
Crop-l-'allow  at  Havre,  Montana 


Available 
Soil  Water 


Probability 


Cumulative  Probability 


(inches) 

0.0  -  1.0 

1.0  -  2.0 

2.0  -  3.0 

3.0  -  4.0 

4.0  -  5.0 

5.0  -  6.0 

6.0  -  7.0 

7.0  -  8.0 


.04 
.10 
.22 
.28 
.22 
.10 
.03 
.01 


.04 
.14 
.36 
.64 
.86 
.96 
.99 
1.00 
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higher  than  under  an  alternate  crop- fallow  rotation,  but  also,  the  land  is 
in  crop  during  seasons  in  which  available  soil  water  at  seeding  is  relatively 
high.   This  fact  is  particularly  important  for  saline  seep  control.   A  three- 
year  rotation  with  fallow  followed  by  two  years  of  crop  would  also  result 
in  a  frequency  of  crop  of  two  years  out  of  three.   But  this  fixed  rotation 
would  be  inferior  to  a  decision  rule  given  by  M  =  3  for  purpose  of  saline 
seep  control.   Soil  moisture  would  never  exceed  three  inches  prior  to  a  year 
of  fallow  under  the  rule  given  by  M  =  3.   This  certainly  would  not  be 
true  for  the  fixed  rotation. 

Basic  Trade-off  Equation 

An  application  of  the  basic  trade-off  equation  (BTE)  was  made  at 
Havre  to  see  if  it  yields  results  consistent  with  the  Markov  chain 
model.   The  BTE  is  stated  as : 

(net  return  from  the  seeded  crop  in  relation  to  available  soil  moisture) 
=  (average  net  returns  from  a  crop  after  a  year  of  fallow) 

-  (costs  of  the  fallow) 

-  (average  annual  net  returns  under  the  best  fixed  crop  rotation) 
At  Havre  the  best  fixed  crop  rotation  is  assumed  to  be  alternate 

crop-fallow  (as  opposed  to  continuous  crop),  and  average  available  soil 
water  is  5.59  inches  after  fallow  and  3.12  inches  after  crop.   Referring 
to  the  expected  conditional  yield  information  for  Havre  and  the  costs 
and  returns  information  presented  in  Tables  19  and  20  the  information 
needed  in  the  BTE  can  be  summarized: 

average  yield  after  fallow  =  16. 72  bushels 

variable  cost  of  fallow    =  $4.00 

variable  cost  after  fallow  =  $20.00 
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Using  a  $3.00  price  for  wheat,  the  right  hand  side  of  the  BTE  becomes 
($3.00  X  16.72  -  $20.00)  -  ($4.00)  -  [($3.00  x  16.72  -  $20.00  -  $4.00)/2] 
=  $30.16  -  $4.00  -  $13.08 
=  $13.08 

Application  of  tlie  BTE  implies  that  if  the  net  return  from  tlie  seeded 
crop  in  relation  to  available  soil  moisture  is  equal  to  or  greater  than 
$13.08,  the  crop  should  be  planted.   Using  the  information  concerning 
expected  immediate  returns  in  Table  20  it  can  be  seen  that  if  the  decision 
is  to  be  made  following  crop,  the  available  soil  moisture  value  which  gives 
an  average  return  of  $13.08  lies  between  3.5  and  4.5  inches.   Interpolation 
yields  a  value  of  3.90  inches. 

If  the  decision  is  to  be  made  following  fallow,  the  available  soil 
moisture  value  which  gives  an  average  return  of  $13.08  lies  between  2.5 
and  3.5  inches.   Interpolation  results  in  a  value  of  3.46  inches.   These 
results  agree  closely  with  the  results  of  the  Markov  chain  model  which 
indicates  that  the  optimal  land  use  alternative  is  fallow  if  the  available 
soil  moisture  level  after  crop  is  less  than  4  inches  and  plant  wheat  if 
the  available  soil  moisture  is  equal  to  or  greater  than  4  inches. 

Mandan,  North  Dakota 

The  expected  yield  relationship  for  Mandan  was  estimated  by  using 
the  Havre-Mandan  yield  equation  reported  in  Table  14  together  with  a 
series  of  weather  data  from  Mandan.   The  expected  yield  estimates  and 
their  standard  deviations  are  presented  in  Table  23. 

Probability  distributions  for  availale  soil  moisture  changes  from  one 
season  to  the  next  under  given  land  use  alternatives  were  estimated  from 
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Table  23 


Conditional  Expected  Yields  of  Spring  Wheat 

For  Specified  Levels  of  Available  Soil  Moisture 

at  Mandan,  North  Dakota 


Available 

Conditional 

Standard 

Soil  Moisture 

Expected  Yield 

Deviation 

(inches) 

(bu) 

(bu) 

.5 

8.26 

8.69 

1.0 

9.43 

8.90 

1.5 

10.70 

9.04 

2.0 

12.05 

9.11 

2.5 

13.46 

9.07 

3.0 

14.87 

9.02 

3.5 

16.25 

8.96 

4.0 

17.60 

8.91 

4.5 

18.92 

8.84 

5.0 

20.21 

8.77 

5.5 

21.47 

8.67 

6.0 

22.  71 

8.55 

6.5 

23.90 

8.43 

7.0 

25.06 

8.29 

7.5 

26.19 

8.15 

8.0 

27.27 

8.00 

8.5 

28.32 

7.84 

9.0 

29.33 

7.68 
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recorded  soil  moisture  at  Mandan.  The  jirocedures  used  at  Mnndnn  were 
identical  to  those  used  i.<t  Havre  wliich  have  been  ex|)lained  earlier  in 
the  report.   Tlie  rc^^ression  estimates  for  Mandan  are  summarized  in  Table  24. 

Again,  there  was  no  soil  moisture  series  available  to  relate  the 
effect  of  fallow  on  soil  moisture  measured  after  a  year  cf  fallow.   The 
estimates  for  the  effect  of  fallow  on  soil  moisture  after  a  year  of  crop 
were  used  in  this  instance. 

Costs  and  Returns  Data 

The  cost  data  derived  at  Havre  was  used  for  Mandan  without  adjustment. 
Gross  returns  and  expected  immediate  returns  were  derived  for  a  100%  yield- 
$3.00  wheat  price  and  a  120%  yield-$3.50  wheat  price.  The  returns  data  are 
presented  in  Table  25. 

Economic  Results 

The  results  of  the  Markov  Chain  Model  are  presented  in  Table  26  and 
graphically  in  Figure  2.   The  optimal  value  of  the  decision  variable,  M  , 
is  zero  which  implies  continuous  cropping.   It  is  obvious  that  there  is 
no  economic  loss  from  continuous  cropping  and  fallow  should  not  be  practiced 
for  purposes  of  soil  water  management  even  if  saline  seep  is  not  a  problem. 
It  is  also  clear  that  the  absolute  difference  in  average  net  returns  (which 
do  not  consider  fixed  costs)  between  continuous  crop  and  alternate  crop- 
fallow  would  more  than  compensate  for  tlie  savings  in  fixed  costs  asso- 
ciated with  a  crop-fallow  practice. 

Basic  Trade-off  Equation 

Application  of  the  basic  trade-off  equation  (BTE)  at  Mandan  yields 
results  consistent  with  the  Markov  chain  model. 
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Table  24 


Estimates  of  Parameters  for  Regression  of 
Available  Soil  Water  in  Time  t  on  Available  Soil  Water 
In  Time  t-1  at  Mandan,  North  Dakota 


Estimates 

Land  Use 
Situation 

a 

e 

a 

FC 

0.84 

0.670 

1.602 

CC 

2.54 

0.424 

1.421 

CF 

4.16 

0.418 

1.315 
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Table  26 


Average  Annual  Net  Returns  in  Relation  to  the 
Decision  Rule  Variable,  Mg  for  Two  Specified 
Yield-Price  Situations,  Mandan,  North  Dakota 


Average 

Annual 

Net  Returns 

Decision 

Rule 

Variable, 

Ms 

100%-$3. 

00 

120%-$3.50 

(inches  available 

soil  water) 

(do 

illars) 

0 

36.40 

60.87 

1 

36.22 

60.62 

2 

35.71 

59.81 

3 

34.70 

57.90 

4 

33.23 

54.82 

5 

31.28 

50.91 

6 

28.97 

46.80 

7 

27.00 

43.89 

8 

25.72 

41.42 

9 

24.85 

40.05 
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L20%  yield-$3.50  price 


100%  yield-$3.00  price 


-L. 


_1_ 


i  A  S         h 

Decision  Kiilc  Var  ialj  I  f  ,  M. 
(inches  of  plant  available  water) 


Figure  2.   Average  Annual  Net  Returns  in  Relation  to  the 

Decision  Variable,  M  ,  at  Mandan,  North  Dakota 

s 
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At  Mandan  Che  best  fixed  crop  rotation  is  continuous  spring  wheat 
and  average  available  soil  water  is  6.22  inches  after  fallow  and  4.45 
inches  after  crop.   Referring  to  the  expected  conditional  yield  information 
for  Mandan  and  the  costs  and  returns  information  presented  in  Tables  23  and 
25  the  information  needed  in  the  BTE  can  be  summarized: 

average  yield  after  fallow    =  23.22  bushels 

average  yield  after  crop     =  18.  79  bushels 

variable  cost  of  fallow      =  $4.00 

variable  cost  after  fallow   =  $20.00 

variable  costs  after  a  crop   =  $25.25 

Assuming  a  $3.00  price  for  wheat  the  right  hand  side  of  the  BTE  becomes 

($3.00  X  23.22  -  $20.00)  -  ($4.00)  -  ($3.00  x  18.79  -  $25.25) 
=  ($69.69  -  $20.00)  -  ($4.00)  -  ($56.37  -  $25.25) 
=  $49.69  -  $4.00  -  $31.12 
=  $14.57 

This  implies  that  if  the  net  returns  from  the  seeded  crop  in  relation 
to  available  soil  moisture  is  equal  to  or  greater  than  $14.57  that  wheat 
should  be  planted.  Using  the  information  in  Table  25  it  can  be  determined 
that  if  the  decision  is  to  be  made  following  a  year  of  crop,  the  available 
soil  moisture  vaJue  that  imiiiics  a  return  of  $14.57  falls  below  the  lowest 
possible  soil  moisture  level.  Therefore,  the  BTE  yields  the  same  results 
as  the  Markov  chain  model. 

Sidney,  Montana  and  Dickinson,  North  Dakota 

Weather  series  were  available  at  Sidney,  Montana  and  Dickinson,  North 
Dakota  and  conditional  expected  yields  were  calculated  at  these  locations. 
This  information  is  presented  in  Tabic  27, 
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Table  27 


Condicioiia]  Expected  Yields  of  Spring  Wlieat  for 
Specified  Levels  of  Available  Soil  Moisture  at 
Sidney,  Montana  and  Dickinson,  North  Dakota 


Conditional  Expected  Yield 
Available 
Soil  Moisture  Sidney  Dickinson 


(inches)  (bu)  (bu) 

.5  8.02  10.44 

1.0  9.28  11.94 

1.5  10.61  13.44 

2.0  11.97  14.92 

2.5  13.34  16.38 

3.0  14.73  17.81 

3.5  16.13  19.22 

4.0  17.52  20.59 

4.5  18.88  21.92 

5.0  20.21  23.21 

5.5  21.51  24.46 

6.0  22.72  25.67 

6.5  24.00  26.84 

7.0  25.19  27.96 

7.5  26.35  29.04 

8.0  27.46  30.08 

8.5  28.54  31.08 

9.0  29.58  32.04 

9.5  30.57  32.95 

10.0  31.53  33.83 

10.5  32.45  34.67 

11.0  33.34  35.47 

11.5  34.18  36.24 

12.0  34.99  36.97 

12.5  35.77  37.66 

13.0  36.51  38.33 
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Probability  distributions  for  soil  moisture  changes  from  one  season 
to  the  next  under  alternative  land  use  decisions  were  not  computed  at  these 
locations  because  of  data  limitations.   Therefore  it  was  not  possible  to 
apply  the  Markov  chain  model.   However,  the  BTE  was  used  to  determine 
the  value  of  the  decision  rule  variable  M  .   A  $3.00  price  for  wheat  com- 
bined with  the  cost  estimates  used  at  Havre  and  Mandan,  resulted  in  a 
value  for  M  which  was  less  than  the  lowest  possible  soil  moisture  level 
at  both  locations.   This  implies  that  continuous  cropping  would  be  the 
most  economic  cropping  strategy  at  both  locations  and  that  if  saline 
seep  control  were  an  objective,  fallow  should  be  diligently  avoided. 
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APPENDIX 

The  wheat  yield  response  surface  and  profit  function  for  wheat  and 
nitrogen  prices  at  $3  per  bushel  and  $0.20  per  pound,  respectively,  are 
given  here  in  tabular  form.   In  addition,  the  slopes  of  the  response 
surface  with  respect  to  each  water  and  nitrogen,  while  holding  the  other 
variable  constant,  are  given  in  tabular  form.   (Technically  these  slopes 
are  partial  derivatives  of  the  yield  equation  with  respect  to  each  water 
and  nitrogen.  ) 

The  table  headings  are  self-explanatory  except  for  columns  four  and 
five  where  the  abbreviated  headings  are  as  follows: 

DYDN  =  slope  of  the  yield  response  surface  with  respect  to  nitrogen 
variation  while  holding  water  at  the  fixed  level  in  column 
one  and  nitrogen  at  the  level  given  in  column  two. 
DYDW  =  slope  of  the  yield  response  surface  with  respect  to  water 

variation  while  holding  nitrogen  at  the  fixed  level  in  column 
two  and  water  at  the  level  given  in  column  one. 


A- 2 


opofvjirr^(T'»-«m>cooorvjtrf^  tr«-»rovooof»>of»>of\j 


fnrororommmfu<\j(\jfV)Ojrvj(\)'»->«-e.-t,-«,-j^^^jj-^ 


u 

CT) 

U-l 

)-l 

3 

CO 

4J 

a) 

rH 

c 

e 

o 

ij 

(U 

D. 

o 

r-l 

<n 

U-, 

XI 

tu 

(0 

Qi 

u 

H 

m 

•a 

rH 

X 

^ 

3 

•H 

(U 

ja 

TJ 

•H 

(3 

c 

>- 

H 

(U 

o. 

4J 

C 

a. 

(11 

M 

< 

1-1 

c 

n  ininir)ir)mtr)inirif?if>ir>ir)ioinmir>irir)irrvfr)0>om 

o  oooooooooooooooooooooooo 

o  (Mofin.i^iof»''.-ivflo<v,fv.'»-«oo-*r^{r  ooN^c^o^  o.'roa.i 

>-  r»>>ocrn\0(ricu-*r^(r»-«m.*>orvo(j00»-»omf^oro 

o  ••••*.•••.......•.••••.. 

lu  •••••.••*.•.•.•.•••••••• 

•-•  cvi(vnja)Oja-'(vj.—  .-f»-i'^»-'w->»-i»-<oooo-*mif>iPir 

z 

UI 

C5 

(T)  OOOOOOOOOOOOOOOOOOOOOOOO 

o-  •.•..•••................ 

K  tocir)Oir)C>ir)Oif>oinC'ifioinoir)o      inoir.  o 

•-•  •^--i(Vi(Viroro-*-*iri/NC>cr^rvociccio-.           .-i»-trv) 

7 


Q. 
lit 


< 
3 


o  c.  o  c;  C  C:  C:  C-  o  c.  c    o  o  c  o  o  o  o  o  o  c:  o  o  o 

u)  If  ir  ir   IT  u   IT-  ir  ir  u  1/  If  If  IT  ir  If  ir  ir  If  >£■  >o  nc  ^  >o 


A- 3 


E_,  o>'.^ooor^<rivoooir>omir''rtrn^rn{vi«-HOcnoc'»cinm'.-<(r>oc»cfr'.-i(r,  rv^ru(rr»^'r-*cv)(rai 

KH  ^o•-^^ooooJ^lP>o^o^ou>-*(^Jl^rvlop^o^^O"»-^f\JR:■..-•o^l^■^^oo^o^lnlf^-*■^oc-*<^-*o^o 

u,  •>•••••••••.••••••••••..•..•«..•.....■••. 

2  cxjcv;— <ooo^cf/r^^£lf)^  mai'^c-r-i'^ajfV'Ojcnn",  (»>ro(\jOjCVj..-t.r-tO0>oc.  r->£>if;f»ia»oirr^oo 

~  *-*-»--*-*romrofnrooooomroif)if)inir>imf>ir>inir>intnif)ir)inifiio-<--*-*-#  ■♦^^^tnioin 


d    r')(v*fvjfvj.-i»-i«-ioooco«-<'^«-<oo<rioo«)rvvONom*-*r'>ojfU'^ooo»-»«-«<\jmnooir^>o 

>?     OOOOOOOOOOCOOO«-«»-i»-'COOOOOOOOOOOOOOOOOOOO--i^^ 


3   o^om-*-*'-'f^'^^ir>^ojoovr'0{v,o..  orv'-iinvjDvcir.  •-ivooN«-'-^(r'irornir-*fvjCT>of^H,-i 
Q  yo  00  -^  en  if)  t^  <y  onjr)^m>o>or^'^^r^onjir/r^o'^coir>or~^(roo--»(Vjnjfvj<\jcvi»-«-*oo^ 

>H.,         .............................        .......... 


jtrc>(vjr>^-*ir)inx>>oirifiir)^x'^cr:*{rr»>t^»-t^r^c-f\)^vf)f^oroc(rcx'of'fv»r/irrnir^n.i 
u    •■•••••••«••••*»«•>«■••••■*>••••■••■■>•>• 

^  if)>o>cvo>£)>c>c>rvOvo>CNC>CNC'>crvrvofoo<r(r  ooc»-''^»-*»-''^'»-'"-'»-t»-'--».-<— '^-•.-looo^o 
■^  ..-^^•.-^•.-«.r-t^^..-«^.r-*^»-<.^,-<^.,-f,rf^.rf^^funj{Virurvj<\j(xj(Vicu(U(Vfu(VJfVjfv<^ 

w 
u 
o 

QjOOOOOOOOOOOOOOOCOOOOOOOOCOCOOOOOOOOOOOOOO 

^inoiroinoiroinoinoiro      moiroinoiroir'C'inotnoirsoirc.ircuDoir)       mo 
■^(\j(»)(T)-*-*inir.'<ci>or^r>oooo(r.  •-i«-icv;<\jc»>cr-*^if>ic>t>>or-.rvococo^a'00»-t«-'  «-• 


H  O  C  O  C  C  C  O  C  O  C  C  C>  O  C    O  O  C  O  O  O  C.'  o  c  o  c  o  o  o  o  c   o  o  o  o  o  o  o  c  o  o  o 


A-4 


i_i  (ri»-«ooor^(r>irvoro-*o»Hv£irv(Virr'Ooarrimocir)oc:»-if»>if;r^(r'-ioj-*  io«cr^t^oc.  ocooocoooc- 
faiTiCuainjou^omir.  ir'-*^vooro-«<»>»-.oocoNO*oinr^r^in'-ivOoGocvorvi<<;ocoo>cfu>ooooc. 

^CT\»-«fvi(»>rri-*  ID  W)  ir>  in  V'  IT  *  ^  m  (\.i  -r^  a  oc!  t^  ira  fu^  >oooo<vjm^m  •>or^r^ivr^  rv  r^vcu)-* 
'^m>o>o>0\0>Cvo>o>o>o>OvC<^>CvO>OvOvOir)if)if>if>>ONO>OvOf^rvr^rvr^r^r^iN.rvrvfvfvrvrvrv 


2:(»)rnrororororof"{v>{V(\jnjnjr>:furuococcc(r(ro~u"'fucou^'^oo<r'^oo*.^r^^of^roovor^ 
Q  io*mrU''-<o(T\oO'rv'voif)-»-{»>fv.i.r-<oo«-<ojf»)^^rn(\jo(riooNOir^rU'^ooor^»o^rnnjoo 
g^^^^w^oooooooocooo  oooojfvj  a,' at  •-'»-••-''-'•-•— i--t«-'Oocooooo 


Q-*vOON--<f»)>or^o>'^njr^-*mvO\or^r^rvr^fv<if\jiraooro>oocofvj-*ir)rvoci(riO»-'»-t{vja)<vj 
>^.....    .....................    ............... 

'^aj<vjcumf»>roror)-*-*-*-*-*^^-»^^^^^.H^^<vi<\j(vj<vjr)f»>ror>rop?oo-*-*^-*-*-* 


z 
w 
o 

(liOOOCOOOOOOOOCCOOOOOOOOOOOOOOOOOOOOOOOOOOO 

i-i  in  O  If  I   O  IT)    C'  IT   Ci  IT)    C>  IT)    C'  IT    C'  If)    C'lf)    C>  If )    Ci  If)  IT.   C>  If)  Ci  ID   C)  IT    C'  If;    O  ID   O  If)  O  ID    O  Ifi   O  IT 


a: 

f!^  o  o  o  CO  C'  a  c  o  c  c  c  c  o  c  c  o  c.  o  c-  c  o  c  o  o  o  cj  c  c  c  o  o  o  o  c  o  c.  o  e  o  o 

Soconocororooofofa    ororororororofor   ororoforrr    ocrO'fT   CTOf    trcr    cr  a-  cr  (r  cr  C    rrcrcrcr 


A- 5 


^v^oln-*lO(^c^^orvJ^o-♦oo^co^^ot^^oocl^^C'f^-*o«-^^o^coJ>-^o-*^«.(^.  rvj^vooo(ro«-trufvj 
>ofV!^oo>0'^-*inmeo..-i,-r(r.*fvjv  u^^  o  en  m  -^  >ii  c^  <t-  ■>o  ^-^  ■*  -*o-*if)roo^  »-<Ovoo^or^«-<ru 

r^r^rv>oxoso>oivr^rvooooaoooooaoocaoooo(>ooocooooaoooocoo>o>cr^rvooocaooocn(r)(ri(ri(r> 


o 


o^r-.--ifuif)oc«-<^r^  .-H-^f^ooorv  om>Cr(rvmsochoointvjoomaiOjfVJfvjtvj'.-t'.-«»-t..-f»-f»-too 
run  -*  \0  •^(T\r^  ^  •*  rU'^trii^  vo^cvj^^or^ir^  nj  ootvj  •*  in  r-r^  mroi»-<o>i^  mm  »-io\'r^in(») 
oooofnnjrvj(v.io_'fv.rv.i'^^^.-»..-«^ccooooooooommmnfvio,'a.ojoj^w^^ 


•   •    • 


o     ■*•■*•-*-*      •-•»-«'^(uoj<\j<vjcuf»)r'!r'5c»}n>m(»>por»>f»>ro(»)ro<»/n>  »-««-t«-««-i€vj(\jOj<v)tvicuo 


2     rvj  — .  o^  vO  cr  *  oo  cvi  *  vo  r-  r^  ^.  ir  (t.  cr  ir^  c  ^  «    o  oj  m  r^.  tvj  ^  oo  ir^  r^  iTj  cv  oc  r^)  rv  o  rvj  m  m  a.'  o  r. 
>-     nmfnrotvjfVifuojrvirvirvifV'mmmm  f»)f»>r»)fnro(»>rooorooororo(virvjnjfviAj{\jfnnf»}<»><»>f»5ro 


oocnoaooococoocoooooooccooooooooccooooooo 


u 
o 

H     oir  o  tn 
i-i     o  o  •-•  •-• 

2        •-<»-'   w-l  ^-^ 


IP    O  IT    O  IT   O  IT    O  IT)   O  IT)    O  IT   C>  If  1   d  IT    Ci  If)   O  IT!    C'  IT 

'^•-'(Vitvimoo^^iPir  \o>oi^r-  cx>oco~(r  OO'^*-' 


in  Clin  oin  cm  oin  oin  a 


a: 
w 
H     o  o  c  c  o  o  c  c  o  o  o  o  c  c  c  c  o  o  c  c;  e  o  c  c  o  o  o  o  r.  r;  c^  o  c  o  o  o  c  o  o  o  o 

■^     <^c^(^<^  coccccoeccce  ooccc  c:ococc-.  e^w^^^^w^^^^^^ 


A- 6 


o 
OS 


rufvj'.—  •^0{^^^lr>fo•t-^o^OJ•r-^^rvJ^tu-*ooJooc^^^cr.  (rr^a.ivc^rvir.  r^^vCOO'^oocTi 
oinr^<cnj^^»-«in>cimirocr^f»>io^cr--«a'f»>-*cvivO\cmr^r^r^^ir^n(v;r^ojmonfVjvc 

(rtrcr.  {Tcoor^xO^oj  o■>o•^^o■^lno^cv.l^coc,■.-l0^^r)'Ol^ococoGocr^^clf)m•-^oor^  r^io^.  -*(rim 
<ri  (r(Ts  (T.  (T.  crcr  o^(r.  (riON>cr^rvoooooc'0^crcr'00oooooooooooocrvNorvr^ocoo(r 


2    ooocTiCr.  cr«-<«-r«-ifuaj«-<r-v^  or^f)OvO(X!criir:ruoPin»-«oo^.—  ^^^co(r)rv•.-lr^o^o^oo^o 
Q  .•«^(^^v^(Vlc>«-»^^lf^•^v(^^nooo^o(nl-^(n^o-*•^(^.f^■*^Jc>^^lf>^no•-^-*^ooc•--'oo^o^^  oooin 


ooooocoo— '^-j^^  ^rom 


Q   ■.^•.-*fuoJ(Vlrv.lCV/{VJfv^■.-^ofr^oo^•^f»)^^.  i^cr'^cvi-*tr(>o>cr~rvooococr^rv>ovOir)orovcocofvj 


Q  forororonfoororororo 


airufvjcucv)<vjOjfVjfVifu{Vj(\jfv)(vitvj(vj 


^rooofuirrvoccirr^ifc\jiio(Virrfv.i.-iooir  o^vtoooc  ^^l^r^Jooal^coflt^(T•rvlr>»-^o:tof'0^^c^ 
w     •••••••••••••••••■••••    ••••••••••••••••••• 

t-i  rvrvococococ'ooocoGocf^cvi-*>o<xicr'^rn-*mscr--oco>oo«-i»-<«-'«-t'-i«-''»--'«-<ofvj-*r^cr«-'f\j 
?"   f0f»>f»>rof»>07r»>rororor»)nj<vcvjrufv;rnror»>rorornrnro^  ■*  ^  ■*  ■*  ^-*-*  ^^  ^rvjfvjcvjromn 


2 

W 

o 
o 


0000  0  000OOOOO0OOO000OC00O000C5  000000000000 


M  inr  ir  oir>  nircinoir 


iTi  c>  in  o  ID  c-  if^  c>  IT)  C'  If)  o  in  o  If :  o  in  c-  in  o  if •  o  in 


in  C'  in  o  in 

T^  w-i  (\i  C\j 


g  o  c  o  o  n  e  o  o  o  o  o  o  c  c  o  cv  o  o  c  o  o  o  o  c  c  o  c  c  o  o  o  o  o  c  o  o  c  o  o  c  o 

<  .-^  ,^  ,-1  ,-1  .^  .-f  ^  ,-» .-I .-.  »-4  fv  (V.  <v-  n.-  «v  (V'  fv  (V  a'  CV'  o.  <v  a  <v  fL  n.  o.  o.  (V  a  fK.  O.I  rv  o,'  r^  r^ 


A- 7 


r,  («vc^ojor^f^ovO'^r^fMvO'-^-*oO'.--i^-*mr^irocvCiocir)f^m^cncrr^fr'>oinceinr^^vocvj 

iZ  iv^tv>o..-('.-(00i-«(ri-*-*  •.-.rocv.ivCvomiroop'ifroccr.irfr^  -a-r^triocir^(Vjfv.;f^oc  ^  o  m  <i  ^oo 

o  rv.rt^f^ocvj(ririr<i>Cvcir-*aiOocinr^roo>0'rt>0'r-'Uio^t\jifjr^(r'.^<\jnj<VjfU'^coc>orn 

^  (^ooo»-«'»-<«-<'^'^«-<«-<'^'^'^«-*'^co>or^oooom(r»ooO'.--i.^'r-f'r-tfucu<Vicvjcu(vjfvj'.--t'.-»^ 


§     vOCTifvjuiK'.-H-*r^'»^^f^om^'^ooif>c\j''-«»-«'>--''-«ooo(ri(rcricrococoooor^i^f^rrivt^-** 
sh     ajo^'r^-*«H<nyoro.^<oincoofVj'inr^orn-*'^ooinfvj<ri>ofuo^vornor^'-*«-iooir)(vjorovoo\cu 
Q     moj{v.'rofvj'>-»'^«-''>-'OOOoooo-->«-<inir>-*  -*-*m  ro(*)rvjajnjfV(--i.-«^ooooooo.-i 


Q    •-<'^'^.-«(vj(v.>(\j(vj(vjnj(vj(\jru(viru(u<vj'^*  ■^•^■.-«--4,-e^--..^.^.^..-«,-i^,-t.^.^«-« 


J    >ofv:vr(rc'»-a'»^monrv)*ir^  fvjcr  ^  onrr,  oinococ(T\>cnj>roco0>t^-*aNrriir'>cirfv;orr>o 


o 

a: 


oocooooooooooooc-cooooocooooooooooocoooooo 


>-<    o  in  o  in  o  If)  C'  in  o  itj  c<  if>  o  If/  c>  tn  o  IT' 
z    <nro-*-*inir<c^rv»vococ(rcroo«-i»-' 


IT  c>  If;  cm  o  in  o  in  C'  in  o  in  o  in  o  in  Oin  o  in  o 
•«-H'^ai<Vi<»)r>-*^min'<o»or^r^oC'OC'fro>  OO'^ 


g    o  o  c  o  c  o  o  c  o  o  c  o  o  o  o  c  o  c.  o  e  o  c  o  o  o  o  c  o  o  o  c:  c  o  o  c  c  o  c  o  o  e. 

<      mm   rr  cr   (T   rr    (r   fr,  c-  (Ty  (^    (r    (T  m  rr,  C^   fT'.    (r  4    ^   ^    4  ^  ^    ^^^w*^^^^^^w»w*^^^^^^ 


A-8 


o 


ru-*(r.  rorvonirr^(rO'-»»-'«-fOo>acvD-*'^oc^OvO'^orvoc-*-*o>ocrvj»-«-*»-«fn(^ovc<o 
r-so>coc-*>oruro(roi^of^Lr)'rt't-fi^oc-*-if)0'«--tf^t^<^><virrO'>cro^O'^voir.  <rr^<rr^oc-* 

o^o^oor^on0^^oon^oc^tVl^^or^rvr^^v■olf){»^o^v^^nJO^v^o>o•-^^^,  (rruirr^cro-H 
»-«vOf^ooaoo^(r.  oo•-•»-^•»^fVJ(VJ<\iOJ{\JCVJfU{\J  (vjru(Vi»-t-^vor^oooc(r  oo»-i»-<«-irvjf\jfviojr>m 


^nmo^o(^J<T^lo»^oc■*•-^r^r^lC■^o^o<rlLr■fVi<\.■•^cCT\r^'^coo^>flcncvJu^oc«-<•*r^o^o^c•(^f\Jlf>oc 
""•inifiiT)-*  ■*  ■*mmro(vnjf\,'»-<.^»-ioooc  oo^*--'»-'vf)vcir)ir'^  -*--*-f»>mf<orv.'fv.>.-»-t»-io 


s  ors.rr)0'^'»-'{rir(r<vi<»;rr'«-ir^  «-«-*ir)-*fV)ocrvjir.  ^ir/roir>o.ioo(r.iro-*>0'>f).*»->vO(ri»H^ 
Q  *^cviOfu^i/ir-oci  o-->ojrroo^4'*^-*-poroa)-^ocrirvirim,-«ccvj^in>crvoo(rv<r<r.  oo 

«-•  t   »  '-«  ^  <^  ^  w  .-I  W  ^..-t  ^  ^  ^.1  ^  ^     •  ■   •   •  ^  .1 


w 

e> 

o  oooooooooooooooooooo  ooocoocooocoooooooooc 

02  ••>>*>••*•>••■••••••.•••..••.•.•••...•••. 

^  m   in  o  If)  o  in  o  in  C)  IP  c>  IT;  r>  in  c>  in  c  IT  r>  IT.  n  IP  C'  IT.   in  o  in  C'  ir  o  in  o  in  o  in  o  ir  o  in 
jj.'^    •-•.-iojruc»)rri-*-*mir>o>orvr^oc«icr<r\oo«-'»H     »-i»-<cvjrv,f»;f»>-*^inirvOvcr^»^ 


w  o  o  o  o  o  o  c:  c  o  c-  o  o  r;  c-  o  cr  o  c  o  I  ■/  c  o  f  3  c  o  Ci  c  cj  o  r  >  o  o  o  c  o  o  o  c  c  o  o 

l^  -♦  in  in  If;  ir  If  m  If  ir  IT  If  ir  ir  u   ir  if  u  u  <  if   in  u-  in  ir  if  if  »£  ^c  >c  >r  vc  >c  •>€  ^c  >c  »c  •>€  <f^  >c  vo  vc  >o 


A- 9 


•r^oo^•r-<^^mc^lr>oooru^^lr!lf>-*■^^•^cn^or^7^^{T^•*ooo■*-<-^^oo■r-»^o^coc^f\JOCVIoc 
2fnrofufvJC\j(ufvj^oo«-<»^rucvjfvjroromrotnmroooa!<\j<u»-»oo«-i"-t(U(V!(urncomoof»)roro 


Q  in«-«c\j>o<^rnirvooo^o>ofuoo-*o>ofvjoo^ororvi..^in(rr»>--<ir^rooo^<riir)«--i>otvioorooir> 
>>  ooooo.-t'^cv)  -**  ^r>ro{\.i(vjf\j.rt.»-«coooO'^'^'r-tojir  -♦-*  nfrojfvjfvj»-«'r-io  ooo 


t  » 


QC^l^'0^oLrlrfnoon^c•ccc^ocl^■r-<^^^^v^Cl^r(^(»)U"l^^^■*<v!t^fvoco^c--!m^^fVJortv*l^ 
>-.o\<r(roo»^<oio-*  oO'.-»cv:f»>-»ir)ifiir)ir'jr>in-**  mrv;'r-tC'^<»-fvj--«ooo«-t^-^--i«-iooo 


JfV'frrr.-fv.-^^^fVl^^r^^.-*(^•r^:»rr),_^vO(^.--^o^^fVJln^c^cOl^)^voc^c^(v;^oof>rv^(T-f\J^ofVJ 
^-lCT^(^(^(^ococr^^cl^r^(^•-<^o■*^r^ocl^<^<^ooc^<rocr^»Clnr«(^»-•^^lr)>f^f^oco^<^ooo 


o 

CtiOOOOOOOOOOOOOOCCOOOOOOOOOOCOO  oooooooooooo 

^  c  If)  o  iT  c  in  o  in  IT  o  If)  C'  IT  C'  in  C'  lo  C'  in  c  in  C'  ir  o  in  o  in  ir  C"  in  d  if i  o  in  o  in  C'  m  o  in  o 
-soooctro^  oo«-''^t\;rom-*  -*mir/vO>cr-t^ocioc(rcr  oo»-«»-<cvjrr-t»/-*-*uiir.»c>>«rv»^  oooocr 


(J  o  C'  c  c  c  o  o  o  o  c  c  c:  c  c  c  c  o  o  c  o  e  c  o  o  o  c  o  o  o  c;  o  c  e  o  c  o  c  o  o  c  o 
<«o>o>c>c>c>c>o>cr^r-r^r^f^r^i^ivr^f^t^Kr^r.r--h-f^r^fvcicororororococeeococecoC'orec 


A-10 


o 

Di 


m-o  (xj--<-»-  fuO'T-tio  rufvj  ^  oc\  otf>fviru>x3  rvj  ■.-<r»)oo>oif>>c  ovOin>ooNir>  -♦if>oo*njf»)>0'^(ri 
(»)(V)fu{U'>-<a\0'.-»«Hru(\iojnrommrnro(rif\jojfU'^'^ONOO'^»-»(V(\Jcvimmrororo(\j<\jtxj'^ 


o     <nrooo(vi'0-*<riir)0tr)  0'>o«-i>o(\j'r^<vjr>j>o»-»NOOino>oi'-i\0''-t>0'»-<>o»-tvO»-«vOi»Hrooorooom 
>H    o-^'^oj(\ju~i^^-*r^oncvjcV'^'^oooO'.-t'^cviojrnir/ir^-*roonfV'Cv;'^.r-ioooo»-»»-<rvj 


I   »   • 


Q    .-i<vj»»)iO'Crv>etn^mrnnnjfviroro<nw#if>>crvociocuoo>aoocr-rvrvorvrvrvoc{ro— tfum 

I       •      I  t      I      t       t 


^ 


oofx»*-*fvj-*o*if)*  ^ini^vc-*  oo^^o   ^  rvocvrrvjirrvjrvtrooir  ooj^oor^in^«-'inr^  r^ 


M   cri<roci^vO-*r^cr»-imir>or^«>ccr<rootr{rorf^vcir)n>oocooj-*>orvococ<r(r(r(rocf^>o 


o 


ooooo  ocoo  oo  ooooooooooooooooo  oooo  ooooooooo 


►-<  in  o  IT  c:  IT  IT  c>in  o  in  o  m  c  in  o  in  o  in  o  in  ci  in  r-  in  ir.  c  m  o  in  c-  ir  o  in  o  ir  C'  in  o  m  o  in 


go  o  ooo  o  ooc  c.  coo  o  oooo  cocoeooccocoonooooo  oooo 

<coooooecoc<r«r<rtr(r(r(rcr<r(r<rr'tr{r<r<r<r(r<rciccocooccccooocoo 
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H   cri>c-*(^njif)ooo>oooo>f^-*»--»r^(\jr^'.-i^xOoc'Oininooo«-<oo(T.  ir)ircr.r^  orvoommrvoj 

o 


OS 


i»ioo«->(riif;«-«>oo-*>coo{r.  c^o^fv.l^(Vlao^oo(Vl'^-*»-•l^*VJ^oofVi-*mlf)•*rvJO«^)(VJ^^.-^-*■»-t 


ix,'^o(^(r,  o•<-^«^(VJfVJrucunJrucucunJnJ•^■i^^oooo<^^oo-t-^•.^ru{\JnJa)(VJ(Vl(\J(V!^•.-»ooo^oo 


ooovOfvjocm«-ir^mcr  vorvjnjvOomr^-rHtnfr  oruinoO'^-^f^crojir.  octr:vo*'^ooir)(Vi0".  r^ 
z   o>-*crirooocv>rvtu>0'^inoinoio»^>o«^r^<v»r^»^in{r>rooocvjNOoui<riro'»--»r^r^(r.  ^OvO'-«fvj 
2<Vionir>ir-*-*ro<r'Cvjoj..-'.-^coO'^'^(VfVjf^ro>oio^-*mroojfv)'--<oooO't--.'.H(\irorn^>f) 


irr^ojai-*oo^.HOoru\Ofv,{r-ococ»-«tnoof)rvirr~c>in<\JOc;»->ir;c>NOtr/irvroiri«-ioo(^ 
3ior^-*<»)<\j'»-''^"^'^'»->'»H'-^funj(»i^-or^crotvi»^>o>oifiu^uiinu";Ui\0>or^occr»-«{Vj-*>oooo 


►J  f»"i  00  (v  •-  >fi  a  cr  vo  ^  m  rv.'  cr  {*"•  *  m  (T^  a.'  fv:  c-  iri  or  r^  vf,  a:  in  vc  m  or   o  c^  iT'  or  tr  r^  ^  m  fvj  oo  rvj  cv.'  r^ 
u>«    ••••••••••«••.>    ••••••••••••»••••»••••••. 

^inmcv'iri^  tT'^rrirorvr^oc  oroor^r-  ^oinr^'^cmvo  oc  o(v.I^^lr)ln^o^o^£)^c>om-*fV)»-o^or 

•z 
w 
o 
oooooocooooooooocooococoooooooooooocoooooo 

^oif>ir>C'io  oinoiooin  cir  oir  oin  oinoinir  cun  oin  oio  oiroinc-ir  oin  oin  oinir 
^..-^•.-^ru^o^^»w*^u>lf'>o■«o^~-^^oooo<^o\oo•-^■^AJ^r)^n^^lnl^^£»^cr^l^or'Cco^<^oo•-^•-tcu 


wo  o  c>  cj  c  e  c:  c  o  o  o  c-  c;  o  o  o  c-  e  c;  c  o  o  o  o  o  c  o  c  r   C'  o  o  o  o  c  c  c  c  o  o  c 

<  C)  o  »-  ■^  — '  ^  ^  ^  w  ^  ^  ^  ^  w  -.-•  ^  «-•—•—<  ^  ^  a  a  cv.  cv  <V'  <v.  a  fv^  (V'  cv  rv.'  c  (v  <v>  tv  <v  a  cv:  (\j  n 
■^Aj  «vjcvjrocv>f\jrvj<vfu(V(rvjfv-ro(\jrvjnjfVjrorvjrvj(VJO-'njc\jnj(vrufvjfvjrurofo{\j<uoj(v^ 


A-12 


o 


f^--«tf)ooof\jfvj(v.'«-«osr^rotri-*mnjinr^u"— •— 'mr»)>ocumocf^»-«oco>o>coo><-<oo 
000^  oO'»H'r-t.r-i.^(Vj^.,-»^.^oo<ri(^oor^ocoocr  ooO'r-i..-<.rt^wooocnoooorv 


2   r^>o>oinin-*mmojojCTN(roo^w<Maj(V(riri--»r^  r^io  ir>  ^  oo  ^  c^  ■*  oo  at  o  o  ■*  r^ 
a    >oo-*oooj>oo-*oocurocnvoruoo-*o>o^r^i'^inoonjincrm>covOfuoo»n— •oo^o 
g    inir^mma.nj^oooo--fv:fv.ic'>^  ^voinir^rofr>oj.^-.->coo— •— •(umrf>^ir) 


g    o<^o^ococo^<r.  (TNO«-»f\jpi^inr^<ri'^rn^f»)Ponjainjrnrr>r»5-*iri>crvcroru-*vooo 
Q    oj--*^^— •— «^^fv:<MOjfv>ajcvi{vtcvjf»i<»)f\jfU{V'fXi<V)(vicvjOjfufv.fV]ajfvj(V/mmronr) 


^j   r^-*oo(r.£>«-»ro»-ir^O(r.  ^cc^o^«-(\Jr^f\•l^t^f^^voc  vro^c;r'vcircirofv;r»C'0>«r 
u     ••••••••••••••••••••«•••«•■•■■«»•«.«« 

^    •-•-*oopoajm-**m.*^mr^^coo>CNCcrai-*voooo.^(vnjfv.;aifv.-'Oooh-  -*fvi 

z 
u 
o 
o 
oi   oocoooooooocooooooooooooocccoaooooooo 

;-;   o in  c IT  o  ir.  c  IT)  c  in  c>  IT  Clin  o IT  o  IT  tr   c  in  r-  ir  c:  in  o  in  r.  in  c  in  o  in  o  in  c  in 


w  o  o  c  o  c  o  c  o  o  o  o  o  c  c>  c  c  r :  c  o  r  o  C"/  c-  o  o  c  o  c  c:  o  c  c  o  r .  c,  o  r> 

•-  rue  <\jfvojrut\jrvjOjtv.f\.i  fu«\>(\j(\jfvjojnjrv(\j(\jojfv;aitvirvj(vj<vjfvjcv.irvjfucv.i{Vfvjrv)fVj 


